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Abstract. The results of first principles electronic structure calculations for the metallic 
rutile and the insulating monoclinic Mi phase of vanadium dioxide are presented. In addition, 
the insulating M2 phase is investigated for the first time. The density functional calculations 
allow for a consistent understanding of all three phases. In the rutile phase metallic conduc- 
tivity is carried by metal tig orbitals, which fall into the one-dimensional dy band, and the 
isotropically dispersing bands. Hybridization of both types of bands is weak. In the Mi 
phase splitting of the dy band due to metal-metal dimerization and upshift of the bands 
due to increased p-d overlap lead to an effective separation of both types of bands. Despite 
incomplete opening of the optical band gap due to the shortcomings of the local density ap- 
proximation, the metal-insulator transition can be understood as a Peierls-like instability of 
the dy band in an embedding background of electrons. In the M2 phase, the metal-insulator 
transition arises as a combined embedded Peierls-like and antiferromagnetic instability. The 
results for VO2 fit into the general scenario of an instability of the rutile-type transition-metal 
dioxides at the beginning of the d series towards dimerization or antiferromagnetic ordering 
within the characteristic metal chains. This scenario was successfully applied before to M0O2 
and Nb02- In the d 1 compounds, the du and bands can be completely separated, which 
leads to the observed metal-insulator transitions. 
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1 Introduction 

At a temperature of 340 K and ambient pressure, stoichiometric VO2 undergoes a 
metal-insulator transition, which is accompanied by a structural transition from a 
high-temperature rutile phase to a low-temperature monoclinic phase [jj], || . Since its 
discovery more than fourty years ago, this transition has been attracting considerable 
interest for fundamental reasons, and for possible applications Q|. The latter result 
from the abrupt change in the resistivity over several orders of magnitude, as well 
as from the fact that the transition occurs near room temperature jl], |^]. Despite 
this strong motivation, our understanding of the origin of the phase transition is 
still far from complete. Several models have been proposed ranging from Peierls- 
§, |, g, g, g, §, § to Mott-Hubbard-type g(| [OJ [jJJ scenarios. They stress, 
to a different degree, the role of lattice instabilities, electron-phonon interaction and 
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electron-electron correlations. So far, neither of these approaches has been successful in 
explaining the broad range of phenomena occurring in vanadium dioxide. In particular, 
a complete and generally accepted picture of the physics of this material has not yet 
been developed. More generally, a comprehensive understanding of the rutile-relatcd 
transition-metal dioxides is not available. 

This class of materials has been intensively studied for a long time, due to the 
complexity of the physical properties As is obvious from Tab. [l] this diversity 

Table 1 Properties of transition-metal dioxides with rutile-related structure (from Ref. fllq]). 
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shows up mainly in the 3d series, which comprises, in addition to VO2, the large-gap 
semiconductor Ti02, the half- metallic ferromagnet CrC>2, and the antiferromagnetic 
semiconductor MnC>2 [[l5], [u|. In contrast, except for Nb02, which, like VO2, un- 
dergoes a metal-insulator transition accompanied by a structural transition, the Ad 
and 5d compounds arc neither semiconducting nor magnetic. Nevertheless, there exist 
several members in each group, which display small but characteristic deviations from 
the rutile structure p], [lq, pf. 

Interestingly, the distortions of the rutile structure found at the phase transitions 
of VO2 and NbC>2 are very similar for all the transition-metal dioxides which display 
such deviations, i.e. for V0 2 , Nb0 2 , Mo0 2 , W0 2 , Tc0 2 , and a-Rc0 2 @- The y 
are characterized by two distinct modes, namely a pairing of the metal atoms along 
characteristic chains parallel to the tetragonal c axis, and a lateral zigzag-like dis- 
placement (l8[ |lj, |l9[ |20|| . This leads to the socalled Mi structure, which has a simple 
monoclinic lattice. The unit cell is twice as large as that of the rutile structure. 
NbC>2 is an exception: below the transition temperature this compound has a body- 
centered tetragonal lattice with a unit cell, which comprises 16 rutile cells Due to 
the strong metal-metal bonding within the chains parallel to the rutile c axis, which 
causes the aforementioned pairing, the above six dioxides have the smallest ratios of 
the symmetrized lattice constants cr/clr (Fig. |l|). 
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Fig. 1 Axial ratios of rutile-related transition-metal dioxides. 



The structural similarities of the six compounds mentioned above are in contrast 
to their completely different electronic properties: while the d 1 compounds VO2 and 
NI3O2 undergo a metal-insulator transition, their four d 2 and d 3 neighbours like all 
other Ad and 5d compounds listed in Tab. |l| remain metallic at all temperatures. An 
explanation for this diversity is still lacking. The different trends in the structural 
and electronic properties within the transition-metal dioxides at the beginning of the 
d series led Goodenough to rule out the Jahn- Teller effect as a major source for the 
crystal structure characteristics || . Mechanisms based on Fermi surface nesting could 
be excluded due to the change in d electron occupation. A comprehensive theory, 
explaining both the structural similarities of all the above transition-metal dioxides 
and the diversity in their electronic properties on an equal footing, is still missing. 

The discovery of two additional insulating low-temperature phases of VO2 led to 
further complications. These phases, denoted T and M2, appear on doping or applying 
uniaxial pressure. They are generally regarded as metastable forms of stoichiometric 
VO2 at ambient pressures. Their crystallographic distortions and magnetic properties 
add to the wealth of phenomena observed in vanadium dioxide, and demonstrate even 
more our lack in understanding. 

We have studied the transition-metal dioxides at the beginning of the d series using 
density functional theory based electronic structure calculations p3|. In particular, 
our investigations so far include metallic M0O2 as well as the other metal-insulator 
transition system besides VO2, namely NbC>2 [^3[ ^4). In the present paper we report 
on results of the different phases of VO2. In doing so we complement previous work 
on the rutile and Mi phases by a more detailed analysis of their electronic properties, 
and we present new results for the M2 phase [ p2[ . 

The purpose of this work is twofold. First, by investigating a large class of related 
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compounds within the same theoretical approach, we are able to develop a more com- 
prehensive understanding of the relevant mechanisms. Second, our density functional 
calculations benefit from the fact that they do not depend on adjustable parameters. 
Hence, despite possible limitations arising from the local density approximation, they 
establish a reference, which can be directly compared to experimental data. Finally, 
our results can be used as the starting point for more elaborate treatments, as e.g. the 
LDA+DMFT approach, which only recently has proven very successful in modelling 
the symmetry conserving metal-insulator transition of V2O3 [^5|, |26|. 

The paper is organized as follows: In Sec. an overview of previous results for 
VO2 is presented. In Sec. [| the crystal structures of the different phases are described. 
Sec. H summarizes the gross features of our calculational method. In Sees. ||, [| and 
we discuss the results obtained for the rutile, Mi, and M2 phases, respectively. Finally, 
a summary is given in Sec. |[ 



2 Summary of previous results for VO2 

2. 1 Experimental results 

The important role of the lattice degrees of freedom for the stabilization of the different 
phases of VO2 has been deduced from several experiments, which revealed striking 
differences between the metallic and the insulating state. Together with the symmetry 
change at the metal-insulator transition these findings were taken as evidence for strong 
electron-phonon coupling and a lattice softening at the transition p7[| . Of course, this 
became most obvious from the fact that crystals could even break at the transition 
[|| §. Further more, X-ray measurements revealed large thermal displacements in 
the rutile phase of both stoichiometric and Cr-doped VO2, which exceed the values 
measured in the respective low-temperature phases as well as in the neighbouring 
rutile compounds TiC>2 and Cr02 considerably [^0], |2^]. Rao et al. observed a large 
anisotropy of the thermal expansion in the rutile phase (a|| w 5 - 6 x a_\_) p9| . p0| , 
which causes a strong temperature dependence of the electric field gradient at the 
vanadium site pl[ . In addition, an anisotropic change of the lattice constants at the 
phase transition was reported j$S ||(|: While the rutile c axis increases by w 1% the 
a axis decreases by w 0.6%. At the same time, the volume changes discontinuously 
by -0.044% and the thermal expansion coefficient drops by a factor of three (3^. j30| . 
Ultrasonic microscopy revealed a strong elastic anisotropy in the metallic phase, which 



disappears almost completely in the insulating state 34 . X-ray diffuse scattering 
measurements performed just above the metal-insulator transition by Terauchi and 
Cohen ]35[| indicated a lattice instability at the tetragonal R-point with a wave vector 
parallel to T-R and a polarization vector parallel to the c axis. This result confirmed 
the symmetry analysis by Brews, who showed the rutile-to-monoclinic transition to be 
compatible with a phonon instability at the R-point |36|. Indeed, using a shell model 
to calculate the phonon dispersion of several rutile- type compounds, Gervais and Kress 
found a phonon softening at the tetragonal R-point for VO2 p^j. Furthermore, the 
analysis of the eigenvectors of the R-point soft mode resulted in a displacement pattern 
of the metal atoms equivalent to that of the monoclinic phase |37|. Finally, an oxygen 
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isotope effect was observed for the transition temperature [p8[ . 

The electronic structure of metallic VO2 has been probed by optical and photo- 
emission experiments. According to optical measurements by Verleur et al. the lowest 
unfilled V 3d levels lie about 2.5 eV above the top of the O 2p bands |39|. This 
result was supported by photocmission measurements by Powell et al. [40|. UPS 
and XPS experiments revealed an w 8.5 eV wide occupied band falling into the low 
binding, approximately 1.5 cV wide V 3d band and an w 6 eV wide group of O 2p 
bands at higher binding energies H| HI H, g, f||. According to O K-edgc XAS 
experiments, the unoccupied V 3d bands extend from Ep to 1.7eV and from 2.2 to 
5.2 eV H|, [0], |7|, [f|, H|, H| . The nearly isotropic electrical conductivity, which is only 
twice as large parallel to the rutile c axis as compared to its in-plane value jl9, 50 , has 
been attributed to the existence of at least two different bands near the Fermi energy 
and a rather complex Fermi surface J^J. The same conclusion was drawn from optical, 
Seebeck coefficient and magnetization measurements p9| , 51, For the insulating 



phase, photocmission studies revealed a sharpening and downshift of the occupied V 3d 
bands |^6|, |4l], [l2], |44|, [l5| relative to the metallic phase. XAS measurements showed 
a splitting of the lowest unoccupied 3d state [^7], Q. Finally, the vanadium atom 
pairing is reflected by the magnetic susceptibility, which shows Curie- Weiss behaviour 
above T c and a nearly constant van Vlcck-likc contribution below T c EC], . 



2.2 Theoretical work 

Theoretical approaches to the insulating ground state of VO2 as well as to the metal- 
insulator transition range from phenomenological models to first principles studies. 
In a survey over oxides with octahedrally coordinated transition metals, Goodenough 
pointed out the important role of direct cation-cation interactions, which add to and 
often supersede the cation-anion-cation interactions 0]. If strong enough, the cation- 
cation interactions may give rise to a covalcnt-typc bonding at low temperatures. This 
bonding would go along with a displacement of the metal atom away from the center 
of symmetry of the anion interstice as well as with a spin pairing of the bonding 
electrons. In the case of VO2, the latter would explain the observed Curie- Weiss 
behaviour above T c and the small temperature-independent van Vleck susceptibility 
below T c |5(J. Since those d electrons, which are responsible for the metal-metal 
bonding, do not contribute to metallic conductivity, the phase transition may be from 
the metallic to the semiconducting state ji). Subsequently, Kawakubo proposed that 
the metal-insulator transition might be due to electronic excitations from the bonding 
to the antibonding state |p3[ . Giving arguments for a linear dependence of the energy 
gap on the number of excited carriers Adler and Brooks demonstrated that a system 
with a half-filled conduction band will lower its total energy and enter an insulating 
state by developing either antiferromagnetic order or a crystalline distortion |5| . 

Lateron, Goodenough presented a schematic energy band diagram for both the 
metallic and the insulating phase || fflj. Using electrostatic considerations for the 
effective ionic charges he placed the oxygen 2p levels well below the vanadium 3d 
states. The latter are subject to octahedral crystal field splitting into lower ti g and 
higher e g levels. Furthermore, the t2 g states, which arc located near the Fermi energy, 
are split into a d\\ state, which is directed along the rutile c axis, and the remaining tt* 
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Fig. 2 Schematic band diagrams for VO2. 



states, see Fig. ^J. In the monoclinic phase, metal-metal pairing within the vanadium 
chains parallel to the rutile c axis causes splitting of the dy band into hlled bonding 
and empty antibonding states. In contrast, the tt* bands move to higher energies 
due to the antiferroelectric zigzag-type displacement of the vanadium atoms. As a 
result, a band gap opens between the bonding dy band and the other ti g bands. 
According to Shin et al. the dy band splitting amounts to ~ 2.5 cV, while the tt* 
bands raise by fts 0.5 eV Nevertheless, in proposing this scheme Goodenough 

questioned the covalent-type bonding arising from metal-metal pairing as the major 
source of the phase transition. Instead, he suggested that the metal- insulator transition 
results predominantly from the increased p-d overlap coming with the antiferroelectric 
distortion of the VOg octahedra (6[ 0. It would lower and raise the bonding and 
antibonding tt and tt* levels. The dy states would be split as a consequence of the 
symmetry change. Eventually this would open the optical band gap. Using the general 
band scheme sketched in Fig. ^, Fleam performed numerical calculations for one- 
dimensional vanadium chains and demonstrated the validity of the scenario outlined 
by Goodenough Q. 

A mechanism differing from the molecular orbital picture of Goodenough was pro- 
posed by Zylbersztejn and Mott. These authors attributed the coupled metal-insulator 
and structural transition to the presence of strong electron-electron correlations espe- 
cially in the dy band rather than to electron-lattice interaction pj. However, as Zyl- 
bersztejn and Mott argue, these correlations are efficiently screened by the tt* bands 
in the metallic phase. In the insulating phase, the screening of the dy electrons is di- 
minished since the tt* bands experience energetical upshift due to the antiferroelectric 
displacement of the V atoms. As a consequence, the narrow dy bands at the Fermi 
energy are susceptible to strong Coulomb correlations and undergo a Mott transition. 
This opens the optical band gap. 

Lateron, Paquet and Leroux-Hugon presented a band model for VO2, which aimed 
at incorporating both the electron-electron interactions and the electron-lattice cou- 
pling on an equal footing |T^] . To do so, the authors started from a tight-binding 
representation for the dy and 7r* bands. In addition, they included Hubbard-type 
interactions for both bands as well as a contribution due to the phonons. Using ex- 
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pcrimcntally derived parameters Paquet and Leroux-Hugon concluded that the metal- 
insulator transition is primarily driven by electron-electron correlations. However, 
these correlations would be strongly renormalizcd by the lattice distortion and the 
electrostatic interaction between the d\\ and the tt* electrons. 

First principles studies for rutile VO2 range from calculations using semiempirical 
or overlapping free atom potentials |55|, |5(| to self-consistent state-of-the-art inves- 
tigations [^7[ Q H [)| ^9). Work on the insulating and monoclinic (Mi) phase 
includes a semiempirical calculation p0| , the molecular dynamics study by Wentzcov- 
itch et al. |8|, ||, state-of-the-art investigations using the LDA and LDA+U |2^. plf . 
and a recent model GW calculation |32| . Most of these calculations confirmed the gen- 
eral molecular-orbital scheme sketched by Goodenough. They found the same order of 
bands for the metallic phase and related the change in electronic structure occurring at 
the phase transition to the simultaneous symmetry-lowering distortions of the crystal 
structure. From a maximum in the generalized susceptibility occurring at the tetrag- 
onal R-point, which could be traced back exclusively to the d\\ bands, Gupta et al. 
concluded that the phase transition could be due to formation of a charge-density wave 
accompanied by a periodic lattice distortion leading to the monoclinic Mi structure 
[pri) . Support for the predominant influence of the lattice degrees of freedom on the 
transition came from first principles molecular dynamics calculations by Wcntzcovitch 
et al. j^]. Their variable cell shape approach allowed for simultaneously relaxing the 
atomic positions and the primitive translations. As a result, starting from different in- 
termediate structures these authors obtained the monoclinic Mi structure as the most 
stable one. In addition, a slightly mctastable rutile structure was found, which was 
about 54 mcV per f.u. higher in energy as compared to the Mi structure. Calculated 
crystal structure parameters for both cells were in good agreement with the experi- 
mental data (see Sec. ||). However, the band structure obtained for monoclinic VO2 
contrasted the experimental finding of an 0.6 eV wide insulating gap |39[ |63|, |46[ pT| . 
Instead, semimetallic behaviour with a band overlap of 0.04 eV was found JSfl . This 
result was attributed to the typical failure of the local density approximation to cor- 
rectly reproduce measured optical band gaps. The situation is thus not unlike that 
in small band gap semiconductors as Ge, where use of the local density approxima- 
tion likewise leads to a metallic ground state. Wentzcovitch et al. suspected that, 
by further strengthening the V-V bonds, the gap should eventually open and, hence, 
VO2 could be regarded as a band insulator. This was interpreted as assertion of the 
energy band scheme by Goodenough || |?J . More recent studies using the local density 
approximation confirmed the results presented by Wcntzcovitch et al. in that they 
also found semimetallic rather than semiconducting behaviour |^2| , |59[| . Nevertheless, 
a recent LDA+U approach and a quasiparticle band structure calculation as based on 
a model GW scheme produced the correct optical band gap |6l], |6|] . 

Despite the huge amount of theoretical investigations there still remain funda- 
mental questions. First of all, none of the previous studies has presented a detailed 
analysis of the band structure or the density of states in terms of the orbitals discussed 
by Goodenough. For this reason, the molecular orbital considerations as well as the 
more recent model approaches lack a profound justification. In addition, most of the 
previous studies, including those of Goodenough as well as Zylbcrsztejn and Mott, 
emphasize the role of the tt* bands for the transition. Yet, as will be discussed in 
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more detail in Sec. |3j, the role of the zigzag-like displacement of the metal atoms for 
the metal-insulator transition is still unclear. This is related to the different response 
of the surrounding oxygen atoms. In VO2, the oxygen atoms stay essentially at their 
original positions and the vanadium atom shifts lead to an an increased metal-oxygen 
bonding. This causes an energetical upshift of the tt* bands. In contrast, in M0O2 the 
metal atom displacements are to a large part compensated by the accompanying shift 
of the oxygen octahedron. As a consequence, the zigzag-like displacement of the metal 
atoms does not lead to an increased metal-oxygen bonding. As a consequence, it does 
not generally act as an antiferroelectric mode. Hence, the zigzag-like displacement 
cannot be regarded as a common origin for the destabilization of the rutile structure. 



2.3 The M 2 phase 



The understanding of VO2 and the metal- insulator transition is even more complicated 
due to the presence of two additional phases of this material. They appear on applica- 
tion of uniaxial stress or on doping of VO2 with small amounts of Cr, Fe, Al, or Ga of 
the order of few percent 20, pi] The phase diagram of Cr doped VO2 as given by 
Pouget and Launois |n] is shown in Fig. ||. This phase diagram differs slightly from 
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Fig. 3 Phase diagram of Vi_ x Cr x 02 according to Pouget and Launois The distortion 
patterns of the metallic chains are given in the insets (open circles in the T phase refer to 
the positions of the M2 phase) . 

that given by other authors, which fact has been attributed to differences in sample 
preparation and stoichiometry p0|, |6^ | . Lateron, a similar phase diagram was reported 
for Vi_ x Al x 02 by Ghedira et al. |65|, |66|. According to Fig. |[ Cr-dopcd VO2 enters, 
on cooling, the monoclinic M2 phase in a first order metal- insulator transition p4| . On 
further lowering of the temperature, transition to the triclinic T phase occurs, which, 
for chromium concentrations smaller than 0.015, is followed by a first order transition 
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to the monoclinic Mi phase p0[] . According to X-ray measurements by Marezio et al., 
each of the three low-temperature phases shows a particular distortion pattern of the 



atomic arrangement especially of the characteristic vanadium chains 1 20 . Whereas in 



the Mi phase both the metal-metal pairing and the zigzag-type lateral displacement 
are observed on each chain, in the M2 phase only half of the chains dimerize and the 
zigzag- type deviations are reserved to the other half of the chains. Finally, the T phase 
is intermediate, as those chains, which in the M2 phase dimerize, gradually start to 
tilt, whereas the zigzag chains start to dimerize until distortions in both chains are 
identical and eventually the Mi phase is reached. Hence, the T and Mi phases can 
be regarded as superpositions of two M2-type displacements with unequal and equal 
weights, respectively The role of the latter displacement pattern as the "fundamen- 
tal" one may be understood from the symmetryconsidcrations presented by Pouget 
et al. as well as by Paquet and Leroux-Hugon |64|, [jjl According to these ideas 
the pairing on one chain is intimately connected with the zigzag-type displacement on 
the neighbouring chains: If the vanadium atoms pair 011 one chain, the atoms 011 the 
neighbouring chains, which from the outset are vertically displaced by half the rutilc 
c axis, move off the chain axis towards the apical oxygen atom whose V neighbours 
of the first sublattice have pulled apart. For this reason, the zigzag-like displacement 
may be viewed as being coupled to the metal-metal pairing along the rutile c axis. 

V 51 -NMR measurements by Pouget et al. confirmed the X-ray experiments per- 
formed by Marezio et al.. They revealed a positive Knight shift for the Mi phase of 
Vi_ x Cr x 02, which is similar to that of pure VO2. In contrast, both a positive and 
a negative Knight shift were observed in the M2 phase indicating two different types 
of vanadium atoms J64|. Whereas the positive Knight shift again points to paired V 
sites the negative Knight shift gives evidence for V 4+ atoms with localized 3d elec- 
trons. Finally, as D'Hacnens et al. found from ESR data, the chromium atoms enter 
as substitutional Cr 3+ ions on the zigzag-type chains J67], |(|| . 

Magnetic susceptibility curves taken by Pouget et al. for Vi_ x Cr x 02 showed a 
striking similarity to corresponding data for Fe-doped VO2 as measured by Kosugc 
fl50| . They could be decomposed into a van Vleck contribution for the paired vana- 
dium chains, a Curie- Weiss term representing the Cr 3+ spins and an almost constant 
residual term. Using a model by Bonner and Fisher ]69[| , Pouget et al. were able to 
attribute the latter contribution to noninteracting spin-i linear Heisenberg chains, 
formed by the equispaced vanadium atoms, with an exchange constant of the order 
of room temperature p4j . The M2 to T transition was thus viewed as a bonding or 
dimcrization transition of a linear Heisenberg chain. This interpretation was confirmed 
by the observed decrease of the magnetic susceptibility at the transition, which is ac- 
companied by an increase of the electrical conductivity by a factor of two |66| . From 
the existence of well localized 3d electrons at the V 4+ sites, Pouget et al. concluded 
that the results should be interpreted in terms of a Mott-Hubbard-type picture. In 
doing so, they attributed the insulating gap to electronic correlations |pl| , 

Lateron, Pouget et al. demonstrated that both the M 2 and T phase could be like- 
wise stabilized by applying uniaxial stress along either the [110] or the [110] direction 
in stoichiometric VO2 samples frcfl . Such uniaxial stress can be viewed as suppress- 
ing the zigzag-type displacements on those chains, which have their tilting along the 
direction of the applied stress, while the respective other chains are not affected. How- 
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ever, from the above symmetry considerations it follows that, as a consequence of the 
reduced tilt on one half of the chains, the pairing on the other chains will be also 
reduced. Additional NMR and EPR data revealed that the local symmetries as well 
as the magnetic properties of the two vanadium sites are identical in pure VO2 under 



uniaxial stress and in Vi_ x Cr x 02 70 



As Pouget et al. pointed out, the critical uniaxial stress for appearance of the 
M 2 phase is so small that the free energies of the Mi and M 2 phases in pure VO2 
must be extremely close at temperatures just below the metal-insulator transition 
[(70). The same conclusion was drawn by Pouget et al. for doped VO2. These authors 
pointed to the fact that, although the M 2 and T phases involve substantial changes 
in the vanadium atom positions, they could be stabilized by impurity concentrations 
as low as 0.2% fii) . According to Pouget et al. the M2 and the T phase thus must 
be interpreted as alternative phases of pure VO2 whose free energies are only slightly 
higher than that of the Mi phase of the pure material [Q. As a consequence, the M2 
phase was regarded as a metastable modification of the Mi phase. In contrast, the T 
phase appears as a transitional state, which displays characteristics of both monoclinic 
phases. 



The discovery of the M2 phase with its two different types of vanadium chains was 
difficult to reconcile with all those theoretical approaches, which had explained the 
insulating ground state as originating from crystal structure distortions with a pre- 
dominant influence of cither the pairing or the zigzag-type displacement. Furthermore, 
from the presence of localized id electrons on the zigzag chains and the antifcrromag- 
netic ordering of the local moments within non-interacting linear Hciscnberg chains 
many authors concluded that band theory were unable to correctly describe the physics 
of the M2 and, hence, also of the Mi phase. Instead, electronic correlations were re- 
garded as essential and localized models called for O, Yet, an explanation of 
the physics of VO2 purely in terms of a Mott-Hubbard-type picture might be difficult 
in view of the complex structural changes at the metal-insulator transitions, which lead 
to the same crystal structures as those of the neighbouring metallic transition-metal 
dioxides. 

Although, as Rice et al. have pointed out again only recently p3| , discussion of 
the monoclinic M2 phase of VO2 is essential for a comprehensive understanding of 
the physics underlying the metal-insulator transition, state-of-the-art first principles 
calculations for this phase did not seem to exist. In order to close this gap we included 
an investigation of the monoclinic M2 phase in the present study p2] [. Again, our aim 
is to provide an as broad as possible picture of the physics of the rutile-type transition 
metal oxides using a well defined first principles approach. 

To be specific, since the low-temperature Mi and M2 phases involve characteristic 
but still relative small deviations from the rutile phase we follow the lines of our work 
on M0O2 and NbC>2 and start out investigating the rutile structure [^3[ . The results 
will then serve as a reference for the subsequent investigation of the low-temperature 
phases. 
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3 Crystal structures and local coordinate systems 

3.1 Rutile structure 

The rutile structure of metallic VO2 is based on a simple tetragonal lattice with space 
group PA 2 /mnm {D\\, No. 136) @. The metal atoms are located at the Wyckoff 
positions (2a): (0,0,0), (i, |, 4) and the oxygen atoms occupy the positions (4f): 
±(u,w,0),±(i + u, i — u,~^). The rutile structure is displayed in Fig. |J. According 




Fig. 4 The rutile structure. Large and small spheres denote metal and ligand atoms, 
respectively. 

to McWhan et al. the lattice constants and the internal oxygen parameter are an 
= 4.5546 A, c R = 2.8514 A and u = 0.3001 @. These numbers will be used in the 
calculations. While Westman reported slightly different lattice parameters for metallic 
VO2 |7l| , the numbers by McWhan et al. arc in very close agreement with the values 
given by Marezio et al. for Vo.976Cr .o2402 (u = 0.3004) as well as with those by 
Ghedira et al. for V0 2 doped with 1.5% of Al (a R = 4.5546 A, c R = 2.8528 A and u = 
0.3001) pp| , ^6|. Finally, we mention the parameters resulting from the first principles 
molecular dynamics calculations by Wentzcovitch et al, a R = 4.58 A, cr = 2.794 A 
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and u = 0.3000, which differ only slightly from the experimental values [||. 

As has been already outlined in our previous work on M0O2 p3[ , the rutile struc- 
ture can be alternatively visualized in terms of a body-centered tetragonal lattice 
formed by the metal atoms, where each metal atom is surrounded by an oxygen octa- 
hedron. Octahcdra centered at the corners and the center of the cell are rotated by 
90° about the tetragonal c axis relative to each other. As a consequence, the lattice 
translational symmetry reduces to simple tetragonal. Thus, the unit cell contains two 
formula units. Octahedra neighbouring along the rutile c axis, share edges whereas 
the resulting octahedral chains are interlinked via corners. Each octahedron has or- 
thorhombic symmetry. Yet, the deviations from tetragonal and even cubic geometry 
are rather small for most compounds and still allow for a discussion in terms of the 
latter. There exist two different metal-oxygen distances. The apical distance is be- 
tween metal and oxygen atoms having the same z value. The equatorial distance is 
between the metal atom and the four neighbouring ligand atoms with z = z m etai ±1/2 
©HQ- 

As we will realize in the following sections it is useful to discuss the electronic 
structure of rutile-type compounds in terms of local coordinate systems centered at 
each metal site. These coordinate systems have been already indicated in Fig. |J. Note 
that due to the different orientation of octahcdra centered at the corner and in the 
center of the rutile cell, the local z axes point alternately along the [110] and [110] 
direction. In contrast to the usual adjustment of the x and y axes parallel to the 
metal-ligand bonds we have rotated these axes by 45° about the local z axes such that 
they are parallel and perpendicular, respectively, to the rutile c axis. 

In Fig. ||, we show the angular parts of the metal d orbitals relative to the local 
reference frame of the central metal atom. With the previous choice of local coordinate 
systems the e g states resulting from the cubic part of the crystal field splitting of the 
metal d orbitals comprise the d 3z 2_ r 2 and d xy orbitals. In contrast, the ti g states are 
built from the d x 2_ v 2, d xz , and d yz orbitals. While the d x 2_ y 2 orbitals point along 
the rutile c and the local y axes, i.e. towards the edges of the basal plane of the 
octahedron, the d xz and d yz orbitals are directed towards the faces. In particular, 
the d yz states point along the (100) directions. As a consequence, the d x 2_ y 2 and 
d xz orbitals mediate a- and 7r-type overlap, respectively, between metal sites within 
the vertical chains formed by the octahedra. In contrast, the d yz orbitals have an 
albeit smaller er-type overlap with their counterparts at metal sites translated by the 
vectors (1,0,0). This is due to the above 45° rotation of the local coordinate system. 
This additional rotation interchanges the d x 2_ y 2 and d xy orbitals and adjusts the d yz 
orbitals parallel to the axes of the rutile basal plane. The overlap of both the d x 2_ y 2 
and d yz orbitals with orbitals of the same kind at neighbouring atoms thus connects 
atoms, which are separated by lattice vectors of the simple tetragonal lattice. In 
contrast, coupling between metal atoms, which are located at the corner and in the 
center of the cell, is mediated by the d xz orbitals. These orbitals point to the voids 
between the metal atoms of the neighbouring octahedral chains where they overlap 
with the d x 2_ y 2 orbitals of these chains. 
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(a) d 3z 2_ r 2 



(b) d x 




(c) d x 2 





(d) d x 



(e) d y 



Fig. 5 Angular parts of the d orbitals. 
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3.2 The Mi structure 



The monoclinic Mi structure of stoichiometric VO2 is characterized by a simple mon- 
oclinic lattice with space group P2i/c (C| fe , No. 14) (l§ According to Ander sson 
the lattice constants and the monoclinic angle are au x = 5.743 A, bu x = 4.517 A, Cm 1 
= 5.375 A, and /3mi = 122.61°, respectively 113]. Lateron, Longo and Kierkegaard re- 
ported values a Ml = 5.7517 A, b Ml = 4.5378X c Ml = 5.3825 A, and Ml = 122.646°, 



respectively |19 , which will be used in the subsequent calculations. The crystal struc- 
ture is displayed in Fig. ^. The unit cell comprises four formula units. The metal 




Fig. 6 Monoclinic Mi structure of VO2. 

atoms as well as the two different types of oxygen atoms occupy the general Wyckoff 
position (4e): ±(x,y,z), ±(x, 5 — y, \ + z). The latter are referred to the standard 
choice 

° \ 

c Ml sm(3 Ml (1) 

-CMi cos/3 Ml / 



a Mi,l = 




, a Ml ,3 — 
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for the primitive translations. The atomic positions given by Andcrsson as well as by 
Longo and Kierkegaard are listed in Tables || and ||. Note that all positions have been 

Table 2 Crystal structure parameters of M1-VO2 as given by Andersson (from Ref. Jist). 



Atom 


Wyckoff positions 


parameters 


x y z 


V 


(4e) 


0.242 0.975 0.025 


Or 


(4e) 


0.10 0.21 0.20 


o 2 


(4e) 


0.39 0.69 0.29 



Table 3 Crystal structure parameters of M1-VO2 as given by Longo and Kierkegaard (from 
Ref. 



Atom 


Wyckoff positions 




parameters 








X 


V 


z 


V 


(4e) 


0.23947 


0.97894 


0.02646 


Oi 


(4e) 


0.10616 


0.21185 


0.20859 


2 


(4e) 


0.40051 


0.70258 


0.29884 



shifted by half of the rutile c axis, i.e. quarter of the monoclinic a axis as compared 
to the high-temperature data. 

Again we complement the measured crystal structure parameters by those resulting 
from the variable cell shape molecular dynamics calculations by Wentzcovitch et al. 
^ |J. These authors reported lattice constants = 5.629 A, b^i = 4.657 A, cm 1 
= 5.375 A, and a monoclinic angle (3m, = 122.56°, respectively [||. The calculated 
atomic parameters are listed in Table |j. As for the rutile structure the agreement 



Table 4 Crystal structure parameters of M1-VO2 as given by Wentzcovitch et al. (from Ref. 
§)■ 



Atom 


Wyckoff positions 


parameters 


x y z 


V 


(4e) 


0.233 0.976 0.021 




(4e) 


0.118 0.212 0.228 


o 2 


(4e) 


0.399 0.685 0.293 
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with the measured values is very good. With only few small exceptions all calculated 
numbers are within the range spanned by the two sets of experimental data, which by 
themselves show some variation. 

The close relationship of the monoclinic Mi structure to the high-temperature 
rutile structure becomes visible from a comparison of Fig. ^| to Fig. ^. Furthermore, 
comparing the lattice parameters of the rutile and Mi structures we realize that the 
primitive translations of the latter can be approximately written as 






ajvfi,! ~ I , a Ml ,2 ~ , a Ml ,3 « \ a R . (2) 



This representation emphasizes the prescription of the monoclinic cell in terms of a 
rutile cell doubled along the rutile c axis. 

Nevertheless, as for M0O2 and NbC>2 there exist three distinct differences between 
both structures: i) In Fig. ^ we observe the characteristic metal-metal pairing along 
the rutile c axis. It gives rise to two different metal- metal distances of 2.619 and 
3.164 A, which contrast the average value of 2.851 A assumed in the rutile structure, 
ii) We witness the zigzag-like in-plane displacements of the vanadium atoms parallel 
to the local z axis, hence, parallel to the diagonal of the rutile basal plane. These 
shifts likewise alternate along both the tetragonal a and c axis. As a consequence of 
the zigzag-like distortion two different apical vanadium oxygen bond lengths of 1.77 
and 2.01 A arise. In contrast, the metal atom pairing leads to two short and two long 
equatorial V-0 bond lengths of 1.86, 1.89, 2.03, and 2.06 A. iii) There is a lattice 
strain which causes deviations of the ratios 

e JLph± = 0.9988 , ' 2cMiC ° S/3mi = 1.0096 
from unity. 

Despite the many similarities of the monoclinic structures of M0O2 and VO2 a 
more detailed analysis of the crystal structure data reveals striking differences be- 
tween both compounds. They regard the zigzag- like in- plane displacements of the 
metal atoms and the distortions of the surrounding oxygen octahedra: i) In M0O2 the 
in-plane shift makes an angle of 18° with the local z axis. In contrast, the vanadium 
atoms, in addition to dimcrizing along the rutile c axis, experience in-plane shifts, 
which are parallel to the local z axis, hence, parallel to the (110) direction, ii) Even 
more striking is the different response of the oxygen sublattice to the displacements 
of the metal atoms. In VO2 the oxygen atoms do not follow the metal atom shifts 
and stay almost at their original positions. In contrast, the oxygen atoms in M0O2, 
while hardly moving parallel to the rutile c axis, follow the lateral displacements of the 
molybdenum atoms to a large degree. Hence, in VO2 the lateral displacement of the 
vanadium atoms causes a shift relative to the rather fixed oxygen octahedra such that 
the p-d-\ike metal-oxygen bonding is affected. This is different in M0O2, where the 
molybdenum atom displacements are compensated to a large degree by the simulta- 
neous shift of the surrounding oxygen atoms. In other words, the zigzag-type in-planc 
displacement of the metal atoms does not automatically lead to an (anti)ferroclcctric 
distortion of the whole octahedron, but may still be compensated by a corresponding 
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displacement of the surrounding oxygen cage. Nevertheless, according to the sym- 
metry considerations pointed out e.g. by Paquet and Leroux-Hugon, the zigzag-type 
shift is still coupled to the vertical metal-metal pairing, The previous observation and 
especially the different behaviour of the two dioxides has strong implications for our 
understanding of the whole class of materials. As has been discussed in Sec. |^, the 
lateral zigzag-like displacement of the metal atoms could lead to energetical upshift of 
the it* bands. However, if compensated by accompanying oxygen displacements, this 
antifcrroelectric mode is deactivated and the ir* bands stay at their original position. 
This has indeed been observed in our calculations for M0O2 p3|. We are thus led 
to the conclusion that the antiferroelectric mode can not explain the tendency of the 
transition-metal dioxides to form distorted variants of the rutile structure. For the 
same reason, this mode must be excluded as the driving force for the metal-insulator 
transition, which is intimately connected to the dcstabilization of rutile. These obser- 
vations arc difficult to reconcile with the scenarios sketched by Goodcnough as well as 
Zylbersztcjn and Mott §, @, 



3.3 The M 2 structure 





f a-M 2 \ 








a J/ 2 ,l — 1 








, aM 2 ,3 — 




v S ) 









In the M2 phase, VO2 crystallizes in a centered monoclinic lattice with space group 
C2/m (Cf/j, No. 12) |(3(|. In contrast to the monoclinic Mi structure, this time 
the monoclinic angle distorts the basal plane of the original rutile cell and, hence, the 
setup of the new primitive translations is different from that of the Mi phase. Usually, 
the crystal structure of the M2 phase is specified in terms of the underlying simple 
monoclinic lattice as defined by the primitive translations [^(| 

cm 2 cos (3m 2 \ 
-CA/ 2 sin/?M 2 ■ (3) 
/ 

The lattice constants and monoclinic angle as resulting from the single crystal mea- 
surements on Vo.976Cro.o2402 by Marezio et al. are au 2 = 9.0664 A, &m 2 = 5.7970 A, 
ca/ 2 = 4.5255 A, and (3m 2 = 91.88°, respectively p(J. There exist two different types of 
metal atoms and three types of oxygens, which occupy subsets of the general Wyckoff 
position (8j): ±(x,y,z), ±(x,-y,z), (i, §, 0) ± (x, y, z), and (|, |, 0) ± (x, -y, z). The 
parameters given by Marezio et al., which will be used in the calculations below, are 
listed in Table || Ghedira et al., who performed single crystal measurements on VO2 
doped with 1.5% of aluminum, also found a simple monoclinic lattice with lattice con- 
stants <zm 2 = 9.060 A, bu 2 = 5.800 A, cm 2 = 4.5217A, monoclinic angle (3m 2 — 91.85°, 
and the atomic positions listed in Table |(] |6^, Q . The striking close agreement of 
both data sets supports the interpretation of the M2 phase as a metastablc modifica- 
tion of the Mi phase, which hardly depends on the dopant and the amount of doping. 
As already mentioned, the M2 phase can be likewise stabilized by the application of 
uniaxial pressure along the [110] direction J7(|. Unfortunetely, no crystal structure 
data exist for this situation. 

For the purpose of the present work we prefer using the related centered mon- 
oclinic lattice, which comprises four instead of eight formula units. An alternative 
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Table 5 Crystal structure parameters of the M2 phase of Vo.976Cro.o2402 as given by Marezio 
et al. @. 



Atom 


Wyckoff positions 




parameters 






X 


y z 


Vi 


(4g) 


0.0 


0.7189 0.0 


v 2 


(4i) 


0.2314 


0.0 0.5312 


Oi 


(8j) 


0.1482 


0.2475 0.2942 


2 


(4i) 


0.3969 


0.0 0.2089 


0, 


(4i) 


0.1000 


0.0 0.7987 



Table 6 Crystal structure parameters of the M2 phase of Vo.985 AI0.015O2 as given by Ghedira 
et al. @. 



Atom 


Wyckoff positions 




parameters 






X 


y z 


Vi 


(4g) 


0.0 


0.7189 0.0 


v 2 


(4i) 


0.2312 


0.0 0.5311 




(8j) 


0.1460 


0.2474 0.2865 


o 2 


(4i) 


0.3975 


0.0 0.2284 


3 


(4i) 


0.0980 


0.0 0.7862 



representation of this lattice originates from a redefinition of the first lattice vector as 
1 / a M 2 \ 

ajw a ,i = r (4) 




while slm 2 ,2 arL< i a A/ 2 ,3 are the same as in Eq. (|^). Using this definition and comparing 
the lattice constants of the rutile and the M 2 phases we note the following approximate 
relation for the primitive translations [EOl 




an 

&M2.1 ~ I I , a M2 ,2 ~ I I , a A/2 . 3 re I -a R ] . (5) 
c R J V 2c R ) \ 

They are equivalent to the vectors given in Eq. (||). Indeed, were it not for the lattice 
strain present in both low-temperature structures, which hinders writing the primitive 
translations as integer linear combinations of the rutile lattice vectors, the primitive 
translations of the Mi and M 2 structure were identical (except for a rotation by 90° 
about the Cartesian z axis). This is a consequence of the fact that both monoclinic 
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phases result from the same R-point instability of the rutile phase. This has been 
pointed out by Brews |36| and by Blount (as mentioned by McWhan et al. p7| ). As a 
consequence, we are able to use the same setup of Brillouin zones for both monoclinic 
lattices. This will facilitate our discussion of the band structures below. The close 
relationship between the two monoclinic phases becomes also obvious from Fig. [?], 
where we display the crystal structure of the M2 phase using the unit cell spanned by 




Fig. 7 Monoclinic M2 structure of VO2. 

the lattice vectors (||) and We identify the vanadium atoms Vi and V2 at the 
corners and centers, respectively, of the underlying stacked rutile cells. Oxygen atoms 
of type Oi are located at the apices of the octahedra centered about vanadium atoms 
Vi . They span the equatorial planes of the octahedra centered about vanadium atoms 
V2. In contrast, oxygen atoms of types O2 and O3 are equatorial atoms for Vi and 
apical atoms relative to V2. Since the latter arc shifted along the axis O2-V2-O3 (see 
below) the identity of these two oxygen atoms is lost and two different apical V2-O 
bond lengths evolve. 
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As for the Mi structure, differences between the R and M2 phases arise from the 
distinct local deformations of the VC>6 octahedra. Again, basic entities are i) the metal- 
metal pairing along the rutile c axis and ii) the zigzag-like in-plane displacements of 
the vanadium atoms parallel to the local z axis, hence, parallel to the diagonal of the 
rutile basal plane (see above). As in the Mi structure these shifts alternate along both 
the tetragonal a and c axis. Furthermore, the oxygen atoms hardly move such that 
the antifcrroelectric mode causes a displacement of the vanadium atoms relative to the 
surrounding oxygen octahedron and is not suppressed as in M0O2. 

Yet, in contrast to the Mi phase, the above basic distortion patterns are not equally 
displayed by both chains. As already discussed in Sec. || and as is visible in Fig. [?], only 
half of the chains dimerize (Vi chains) while the other half (V2 chains) experiences 
the zigzag-like displacement without showing any pairing. As a consequence, two 
different V1-V1 bond lengths of 2.538 and 3.259 A appear while there is only one V2- 
V2 bond length of 2.933 A. Within the octahedra centered about Vi the two apical 
V1-O1 distances are 1.868 A. Due to the vanadium-vanadium dimerization within 
these chains two short V1-O2 distances of 1.852A and two long V1-O3 distances of 
2.089 A arise. In contrast, in the octahedra centered about the tilted vanadium atoms 
the four equatorial V2-O1 bond lengths of 1.931 and 1.974A are very similar. At 
the same time, the distances between the central V2 atom and the apical O2 and O3 
disproportionate and assume the values 2.127 and 1.726 A. 

Finally, we mention the lattice strain, which, in addition to the strains present 
in the Mi lattice, distorts the rutile basal plane. The resulting parallelogram would 
be likewise expected to result from uniaxial pressure along the [110] direction. It 
suppresses the zigzag-like displacement on those chains, which have their local z axes 
parallel to the shorter diagonal of the parallelogram. Since, according to the above 
symmetry arguments, the zigzag mode is connected to the dimerization of metal atoms 
on the neighbouring chains, suppression of the in- plane displacement on half of the 
chains may well lead to a suppression of the metal-metal pairing on the other half. 



4 Method of calculation 

The calculations performed in this study arc based on density functional theory (DFT) 
and the local density approximation (LDA) J73|, [74|. As in our previous work on CrC>2 
[H, |6|, M0O2 H!|, and Nb0 2 @] we employ the augmented spherical wave (ASW) 
method in its scalar-relativistic implementation (see Refs. j78[ (79J |o) for more 
recent descriptions). Since the ASW method uses the atomic sphere approximation 
(ASA) pljl , we had to insert so-called empty spheres into the open crystal structures 
of VO2. These empty spheres are pseudo atoms without a nucleus, which are used to 
model the correct shape of the crystal potential in large voids. In order to minimize the 
sphere overlap we have recently developed the sphere geometry optimization (SGO) 
method, which solves the problem of finding optimal empty sphere positions as well 
as radii of all spheres automatically . The routine was applied to all three crystal 
structures under consideration. For the rutile structure, addition of 16 empty spheres 
allowed to keep the linear overlap of any pair of physical spheres below 15% and the 
overlap of any pair of physical and empty spheres below 20%. The positions of the 
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empty spheres are listed in Tabic M. In addition to the empty sphere positions the 



Table 7 Empty sphere positions for the rutile structure. 



Atom 


Wyckoff positions 




parameters 






X 


y z 


Ei 


(4c) 






E2 


(4g) 


0.3238 


-0.3238 


E3 


(16k) 


0.0109 


0.2973 0.2358 



SGO algorithm proposed atomic sphere radii for all spheres, which are listed in Table 
H together with the orbitals used as the basis set for the present calculations. States 



Table 8 Atomic sphere radii and basis set orbitals used for the rutile structure. 



Atom Radius/as Orbitals 



V 


2.230 


4s 


4p 


3d (4f) 





1.831 


2s 


2p 


(3d) 


Ei 


1.572 


Is 


(2P) 




E2 


1.692 


Is 


2p 


(3d) 


E:>, 


0.949 


Is 


(2p) 





given in parentheses were included as tails of the other orbitals (see Refs. [T^, 78, BO) 
for more details on the ASW method). 

Applying the SGO algorithm to the monoclinic Mi structure we were able, by 
inserting 48 empty spheres, to keep the linear overlap of any pair of physical spheres 
below 16%, and the overlap of any pair of physical and empty spheres below 21%. The 
positions of the empty spheres, the atomic sphere radii for all spheres as proposed by 
the SGO algorithm as well as the basis set orbitals are listed in Tables ^ and [To[ 

Finally, for the monoclinic M2 structure, we were able, by inserting 29 empty 
spheres, to keep the linear overlap of any pair of physical spheres below 17%, and the 
overlap of any pair of physical and empty spheres below 22%. The positions of the 
empty spheres are listed in Table |ll| and the corresponding atomic sphere radii as well 
as the basis set orbitals are given in Table [l^. 

Self-consistency was achieved by an efficient algorithm for convergence acceleration 
[ §3j . The Brillouin zone sampling was done using an increased number of k-points 
ranging from 18 to 1800 points, 54 to 6750 points, and 63 to 1088 points, respectively, 
within the irreducible wedges of the three Brillouin zones. This way we were able to 
ensure convergence of our results with respect to the fineness of the k space grid. 

In addition to analyzing the band structure and the (partial) densities of states 
we will also address chemical bonding in terms of the crystal orbital overlap popula- 
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Table 9 Empty sphere positions for the Mi structure. 



Atom 


Wvckoff positions 




parameters 








X 


y 


z 


i 

Ei 


(2d) 


U.o 


U.O 


U.U 






n 1 ^Afi 

U. lO^tU 


n 00^7 

-U.UUO ( 




E3 


(4c) 


U.4UZ0 


U. 14oo 


n 01 no 




(4e) 


U.Uooo 


9/1 ^9 


U. 1DUZ 


E, 


(4e) 


0.3274 


0.2118 


0.4983 


E(i 


(4e) 


0.2726 


-0.1788 


0.4289 


E 7 


(4e) 


0.0141 


-0.0326 


0.2357 


Em 


(4e) 


0.2317 


-0.0336 


0.2812 


Eg 


(4e) 


0.1289 


-0.0016 


-0.2788 


Eio 


(2b) 


0.5 


0.0 


0.0 


En 


(4e) 


0.4375 


0.1054 


0.2706 


El2 


(4e) 


0.1704 


0.2552 


0.0292 


E13 


(4e) 


0.3768 


-0.0322 


-0.3779 



Table 10 Atomic sphere radii and basis set orbitals for the Mi structure. 



Atom Radius/as Orbitals 



V 


2.042 


4s 


4p 


3d (4/) 




1.841 


2s 


2p 


(3d) 


o. 


1.739 


2s 


2p 


(3d) 


Ei 


1.832 


Is 


2p 


(3d) 


E2 


1.753 


Is 


2p 


(3d) 


E:j 


1.772 


Is 


2p 


(3d) 


E 4 


1.734 


Is 


2p 


(3d) 


Eg 


1.236 


Is 


(2p) 




Efj 


1.213 


Is 


(2p) 




E7 


0.978 


Is 


(2p) 




E,x 


0.974 


Is 


(2p) 




Eg 


0.992 


Is 


(2p) 




Eio 


0.841 


Is 


(2p) 




En 


0.794 


Is 


(2p) 




E12 


0.791 


Is 


(2p) 




E13 


0.841 


Is 


(2p) 
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Table 11 Empty sphere positions for the M2 structure. 



Atom Wyckoff positions parameters 









11 




Ei 


(8j) 


0.2491 


0.1251 


-0.0019 


Eo 


(2c) 


0.0 


0.0 


0.5 


E:; 


(4h) 


0.0 


0.3149 


0.5 


E 4 


(4i) 


0.0765 


0.0 


0.1518 


E 5 


(8j) 


0.1597 


0.2525 


-0.3270 


Eo 


(4i) 


0.3736 


0.0 


-0.1461 


E 7 


(8j) 


0.0075 


0.1525 


-0.2867 


E,s 


(4i) 


0.3932 


0.0 


-0.4322 


Eg 


(4i) 


0.2241 


0.0 


0.2109 


E10 


(4i) 


0.4763 


0.0 


-0.3109 


En 


(8j) 


0.1604 


0.2567 


-0.0648 



Table 12 Atomic sphere radii and basis set orbitals for the M2 structure. 



Atom Radius/aB Orbitals 



Vi 


2.213 


4s 


4p 


3d 


(4/) 


V 2 


2.085 


4s 


4p 


3d 


(4/) 


Oi 


1.848 


2s 


2p 


(3d) 




Oo 


1.817 


2s 


2p 


(3d) 




0, 


1.705 


2s 


2 P 


(3d) 




Ei 


1.642 


Is 


■2 P 


(3d) 




Eo 


1.892 


Is 


2 P 


(3d) 




E 3 


1.978 


Is 


2p 


(3d) 




E 4 


1.929 


Is 


2p 


3d 


(4/) 


E, 


1.795 


Is 


2p 


(3d) 




E„ 


1.781 


Is 


2p 


(3d) 




E 7 


1.104 


Is 


(2p) 






E s 


1.195 


Is 


(2p) 






Eg 


1.143 


Is 


(2p) 






E10 
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tion (COOP) as based on the notions introduced by Hoffmann |Q as well as related 
concepts. Among the latter are the crystal orbital Hamiltonian population (COHP) 
as proposed by Dronskowski and Blochl [ ]S7^ and the covalence energy by Bornscn et 
. All these quantities have been implemented in the ASW method [§5| (see 
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also Refs. |l|, and have been validated for a large number of compounds fllrj| . 



5 Results for metallic V0 2 

5.1 Molecular orbital picture 



Our expectations on the electronic structure of VO2 are easily stated within the molec- 
ular orbital picture proposed by Goodenough (see Sec. ^2) || ||, 0. The general situ- 
ation is visualized in Fig. ||. First, hybridization between the oxygen 2p and vanadium 
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Fig. 8 Molecular orbital scheme for rutile VO2. 



3d orbitals will lead to a- and 7r-type overlap. This gives rise to states of a and a* as 
well as 7r and n* character. Since the p-d overlap is stronger for the former these states 
experience the larger bonding-antibonding splitting. While the a and 7r states will be 
filled and primarily of O 2p character, the corresponding antibonding bands will be 
dominated by the V 3d orbitals. Due to the electron count, the latter will be found 
at and above the Fermi energy. In the nearly perfect cubic octahedral surrounding of 
the vanadium atoms, the a* and 7r* states arc of e g and ti g symmetry, respectively. 
We identify these orbitals in Fig. ^ as the d 3z 2_ r 2 /d xy and the d x i_ y i j d xz j d yz or- 
bitals, respectively. Note that the albeit small orthorhombic distortions of the oxygen 
octahedra lift all the degeneracies among the d orbitals. As already pointed out by 
Goodenough, the exact position and width of the d bands is subject not only to the p-d 
hybridization but also strongly influenced by direct metal-metal interactions j^, [| . In 
particular, in the rutile structure such interactions involve the d x 2_ y 2 orbitals, which 
experience strong overlap parallel to the rutile c axis. These orbitals are of bi g sym- 
metry but are usually designated as the d|| bands. The remaining t 2 g orbitals are of e g 
symmetry and usually subsummed under the name n* (not to be mixed with the same 
notation used above for the t2 g states as a whole). While the d xz orbitals experience 
hardly any overlap with orbitals of neighbouring metal atoms the d yz states mediate 
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metal-metal bonding parallel to the in-plane axes of the tetragonal cell. The general 
situation was already sketched in Fig. ||. 

The validity of the just outlined general band scheme for the rutile-type transition- 
metal dioxides has been already confirmed in our previous calculations for hypothetical 
rutile M0O2 as well as metallic NbC>2 [^3[ Furthermore, from the close relationship 
of the crystal structures of all three compounds we expect a very similar electronic 
structure also for VO2. Yet, as compared to the Ad systems the vanadium d bands 
will display smaller bandwidths. 



5.2 Band structure and density of states 

We display in Fig. ^ the electronic states along selected high symmetry lines within 




Fig. 9 Electronic bands of rutile VO2 along selected symmetry lines within the first Brillouin 
zone of the simple tetragonal lattice, Fig. D-OKa). 



the first Brillouin zone of the simple tetragonal lattice, Fig. 10(a). The corresponding 
partial densities of states (DOS) are given in Fig. O. The total density of states at 
the Fermi energy, N(Ep), is 2.63 states/f.u./eV. Not shown are low lying oxygen 2s 
states. 

In Figs. |9] and [ll] we identify four groups of bands. In the energy range from -7.6 
to -2.2 eV we observe 12 bands, which trace back mainly to O 2p states but have a 
non-negligible contribution due to the V 3d states. Bands are most easily counted 
along the direction X-R where they are twofold degenerate. The next two groups, 
which extend from -0.6 to 2.0 eV and from 2.0 to 5.5 eV, contain six and four bands, 
respectively. They originate mainly from V 3d states. Yet, p-d hybridization causes 
additional O 2p contributions in this energy range. From the above molecular-orbital 
point of view we interprete the bands between -7.6 and -2.2 eV as the O 2p-dominatcd 
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(a) (b) 

Fig. 10 First Brillouin zones of the (a) simple tetragonal and (b) simple monoclinic lattices. 
Backfolding of the tetragonal Brillouin zone implies the following transformation of high 
symmetry points: Xt — > Ym, Zm; Mt — > Cm; Zt — » Ym; Rt — » T; At — ► Zm- 




(E-Ep) (eV) 

Fig. 11 Partial densities of states (DOS) of rutile VO2 per formula unit. 
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lower a and higher 7r states. The respective 7r* and a* states are found in the energy 
intervals from -0.6 to 2.0 eV and 2.0 to 5.5 eV. Finally, we observe V 4s states starting 
at 7.0 eV. 

Crystal field splitting expected from the fact that the metal atoms arc located at 
the centers of slightly distorted V06 octahedra is observed in the partial V 3c? DOS 
shown in Fig. U2. There we have included only the single vanadium atom at the corner 
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Fig. 12 Partial V 3d t2 g and e g densities of states (DOS) of rutile VO2. Selection of orbitals 
is relative to the local rotated reference frame. 

of the rutile cell and used the local rotated reference frame introduced in Sec. [Tl]. Fig. 
fl2| clearly reveals the almost perfect energetical separation of the 3d ti g and e g groups 
of bands. The former states appear almost exclusively in the energy range from -0.6 to 
2.0 eV. The e g states dominate the bands between 2.0 and 5.5 eV. The small but finite 
t2g~£g configuration mixing is a measure of octahedral distortions. Contributions of 
the V 3d states to the oxygen bands are slightly larger for the e g states which, forming 
a bonds, experience a larger overlap with the O 2p states. For the same reason, the 
bonding-antibonding splitting is larger for the e g states as compared to the ti g states, 
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which give rise to n bonds. 

In Fig. 12 the differences between the single symmetry components of the t2 g and 
e g orbitals, respectively, reflect the orthorhombic site symmetry. Worth mentioning 
is the pronounced double peak structure of the d yz partial DOS. Since the respective 
orbitals point along the two in-plane a axes of the tetragonal cell these double maxima 
result from metal- metal bonding. Finally, the d x 2_ y 2 partial DOS displays an albeit 
smaller separation into two peaks, which is due to the overlap of these orbitals parallel 
to the rutile c axis. 

Our findings are in perfect accordance with the molecular orbital picture proposed 
by Goodenough j^, J7|, which was sketched in Sec. 5.1. There the d x 2_ y 2 states were 
designated as the d\\ states and the notation 7r* was used for the d xz and d yz states. 
Deviations from the band diagram of Goodenough show up with respect to the widths 
of the different bands. Goodenough and many authors after him assumed the dy 
band to be considerably more narrow than the tt* bands. Yet, our calculation reveals 
similar widths of all three ti g bands. In particular the dy band has a width of ~ 
1.5 eV. Of course, this has strong implications for model approaches building on strong 
correlations within this band. 

While the results for rutile VO2 agree very well with those for hypothetical rutile 
M0O2 |23| differences appear in a more detailed analysis of the single symmetry com- 
ponents of the t2 g orbitals as given in Fig. |l2|. In M0O2, the Mo Ad x 2_ y 2 band shows a 
strong tendency towards bonding-antibonding splitting while the remaining t2 g bands 
give rise to very similar densities of states. In contrast, in VO2 bonding-antibonding 
splitting of the V 3d x 2_ y 2 band is much weaker. At the same time, the 3d xz and 3d yz 
partial DOS show distinct deviations. In particular, the latter reveals some splitting 
into two peaks at about -0.3 and 1.3 eV. These features may be traced back to the 
much larger c axis of VO2 (2.8514A vs. 2.805A in M0O2), which leads to a larger 
separation of the metal atoms along this axis. In contrast, the rutile a axis is reduced 
(4.5546 A vs. 4.856 A in Mo0 2 ). This causes a larger overlap of metal d orbitals within 
the planes and, hence, a visible bonding-antibonding splitting of the local d yz orbitals. 

Whereas our calculations are in good agreement with the calculations by Caruthers 
et al. as well as by Gupta et al. (5f| a closer look reveals small but distinct 
deviations concerning the exact band positions and, consequently, the Fermi surface. 
We attribute these differences to the approximations inherent in the old calculations. 
Good agreement, however, is found in a comparison with the existing state-of-the-art 
calculations by Wentzcovitch et al. || as well as Kurmaev et al. |5{| . Minor differences 
in the band structures may result from the use of the crystal structure data of Ref. 
|f7l| or from Brillouin zone samplings with different numbers of k-points. Finally, we 
observe distinct deviations in the recent calculations by Nikolaev et al. with respect 
to the band positions ^7[. While our calculations give an occupied V 3d bandwidth 
of 0.6 eV, in close agreement with Wentzcovitch et al., Nikolaev et al. report a value 
of 0.77 eV. At the same time, the latter authors locate the upper edge of the O 2p 
bands 1.39 eV below the Fermi level. In our calculation the band maximum is found 
at -2.4 eV in almost perfect agreement with both the optical measurements by Verleur 
et al. and the photoemission data by Powell et al, Shin et al., Bcrmudcz et al. as 
well as Gocring et al.. These authors find the O 2p states 2.4-2.5 eV below the Fermi 
energy |9[ |fj, |l|, ||, || || . The total width of the O 2p band of 5.37 eV reported 
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by Nikolaev et al. is again in agreement with our value of ~ 5.6 cV. Finally, we point 
to the almost equal densities of states at the Fermi energy of all three tig orbitals, 
which confirm the reported isotropic conductivity |49, |3Cj . 



5.3 Chemical bonding 

Chemical bonding is addressed via the covalence energy shown in Fig. [l3| The curves 
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Fig. 13 Partial covalence energies (E cov ) of rutile VO2. 



display the "canonical" behaviour being negative (bonding) and positive (antibond- 
ing) in the low and high energy regions. While the oxygen-oxygen overlap is almost 
negligible we observe considerable metal-metal bonding in the V 3d dominated groups 
of bands. Obviously, the dominating contribution to the total bonding results from 
the V 3d-0 2p bonding. Below Ep, the corresponding E cov curve is negative below 
-2.2 eV and positive only in the small energy range of the occupied V 3d derived bands. 



5.4 Comparison to experiment 

We display in Figs. |lj and [IB] total and partial V 3d and O 2p densities of states folded 
with a Gaussian of 0.25 and 0.5 eV width, respectively, for the occupied and unoccupied 



part of the spectrum. In Fig. 14 we have added UPS spectra as measured by Goering 
et al. |45|j . Good agreement of the calculated and measured curves is found. In 
addition, the calculated total DOS compares very well with the XPS and UPS spectra 
by Blaauw et al., Shin et al. and Bermudez et al. Jl6| [II], [l2). The occupied bandwidth 
of w 8.5 eV deduced from the experiments is close to the calculated 8.0 eV. According 
to the XPS and UPS spectra the valence band is split into a low and high binding 
part of about 1.5 and 6eV width, respectively M, Ell Ea|. Whereas the low binding 
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Fig. 14 Total and partial densities of states (DOS) of rutile VO2 folded with a 0.25 eV wide 
Gaussian (lower set of curves) and UPS spectra (upper curve; note the offset introduced in 
order to distinguish experimental and theoretical results; from Ref. j43L Eq|). 
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Fig. 15 Partial O 2p densities of states (DOS) of rutile VO2 folded with a 0.5 eV wide 
Gaussian (lower set of curves) and XAS O K edge spectra (upper set of curves; note the offset 
introduced in order to distinguish experimental and theoretical results; from Ref. PQ , Q; 
data shifted by 529.5 eV). 



V. Eyert, The metal-insulator transitions of VO2 



31 



region results from a single peak centered at w — l.OeV, which is attributed to the V 3d 
states, the high binding region, resulting mainly from the O 2p states, is dominated by 
a large peak at m — 5.5cV and shows two slight shoulders at -4.0 and -7.5 eV |^3|, ff5[ . 
In contrast, as compared to the UPS data, the XPS spectra display a considerably 
higher intensity at energies below -5.5 eV and a strong decrease in intensity around 
the top of the O 2p band p6l [H| . Since the V 3d states have a larger cross section in 
XPS ]9l| this shift of intensity reveals a considerable V 3d contribution especially to 
the bands below -5.5 eV and, hence, a strong p-d hybridization in this energy region. 
This is well reflected by the calculations. 

In Fig. |l5| we have complemented the theoretical results by soft-X-ray absorption 
spectra as measured by Muller et al. [Bo| The experimental data are shifted by 
529.5 eV. In order to probe the angular dependence, in these experiments the polar- 
ization vector E was oriented either parallel (4> = 90°) or else perpendicular (0 = 0°) 
to the rutile c axis. In the former case, dipole selection rules allow for transitions from 
O Is to the O 2p z state. In contrast, for E perpendicular to the c axis transitions to 
the O 2p x and 2p y states may occur. The curves given by Muller conform with the 
spectra by Abbate et al. who found, for metallic VO2, two peaks at 529.9 and 532.5 eV 
photon energy for E || c jl7). In the spectra given by Muller et al. these peaks arc 
located at 530.3 and 532.7 eV, respectively 0]. However, as Muller has pointed out 
previously, the curves actually can be fitted to four peaks at 530.1, 530.7, 532.7, and 



534.7 eV |90|. In the energy-shifted spectra of Fig. |15| these peaks are located at 0.6, 
1.2, 3.2, and 5.2 eV. 

In Fig. [l5|, the partial DOS correspond to the final states, which are occupied 
in the XAS experiments. Note that the curve marked O 2p Xl p y comprises only the 
mean average of these orbitals since the component perpendicular to E is not seen in 
experiment. For the same reason, the curve marked O 2p actually contains the full 
contribution from the p z orbital but only half of the contributions from the p x and p y 
orbitals. As for the occupied part of the spectrum we observe overall good agreement 
between experiment and theory. This holds for the positions of the peaks, their relative 
intensity and their angular dependence. Still, agreement seems to be less satisfactory 
in the high energy region between 4 and 5 eV. However, the discrepancy is resolved by 
the aforementioned four peak analysis of the experimental spectra. Finally, differences 
with respect to intensities may be attributed to matrix element and core-hole effects, 



which might be important but are not accounted for in our present analysis 92, kJ\ 



5.5 Fermi surface 

Despite the good agreement of our results with previous experimental and calculated 
data we are still seeking for more conclusive arguments for the destabilization of the 
rutile structure of VO2 at low temperatures. Since the metal-insulator transition is 
accompanied by a structural transformation one might be tempted to deduce more 
information about the origin of the transition from investigating the Fermi surface. 
This has been already done by Gupta et al. in the course of their non-selfconsistent 
calculations |5(|. Furthermore, investigation of the Fermi surface of hypothetical rutile 
M0O2 in our previous work revealed strong nesting on flat portions at height q z = 
ij-- This indicated a possible splitting of the d x 2_ y 2, which was indeed observed in 
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monoclinic M0O2 p3| . 

Cuts through the Fermi surface parallel to the {100} and {110} planes as resulting 
from the present calculations are shown in Figs. [l6[ [H], and [l8| They contain the F-X- 
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Fig. 16 {100} cut through the Fermi surface of rutile VO2. 

Z-R, T-M-Z-A, and X-M-R-A planes, which are the vertical borders of the irreducible 
wedge of the first Brillouin zone, Fig. |l(](a). As before, interpretation of the figures 
is greatly facilitated by combining them with the band structure shown in Fig. ||. 
We recognize electron- and hole-like regions near the T- and Z-point, respectively. 
Nevertheless, in contrast to the non-selfconsistent calculations by Gupta et ai, our 
calculations do not allow to identify any Fermi surface nesting in Figs. [l6]to[l^, which 
might be indicative of a destabilization of the rutile structure. This is in contrast to 
the aforementioned findings for hypothetical rutile Mo02- However, since VO2 has one 
d electron less, the differences between both compounds were to be expected. We thus 
conclude, that investigation of the Fermi surface does not allow for an explanation 
of the instability of the rutile structure. In so far, our result resembles the above 
mentioned argument of Goodenough against a Jahn- Teller-type distortion as the major 
source for the crystal structure deformations. 

5.6 Symmetry analysis of the band structure 

In order to find a common origin of the instability of the rutile structure in the d 1 , 
d , and d? transition- metal dioxides we proceed along the same lines as in our work 
on M0O2 and Nb02 and perform a detailed analysis of the electronic states. To this 
end we display in Figs. |l^, |2^, and ^l] the electronic bands in the narrow energy range 
from -1.0 to 2.5 eV about the Fermi energy in the same special representation as in 




V. Eyert, The metal-insulator transitions of VO2 



33 




-0.4 -0.2 0.2 0.4 

k (1/Bohr) 



Fig. 17 {110} cut through the Fermi surface of rutile VO2 



0.4 



0.2 



eg - 

M 



M 



-0.2 



-0.4 



-0.3 -0.2 -0.1 0.1 0.2 0.3 
k (1/Bohr) 



Fig. 18 {100} cut through the Fermi surface of rutile VO2. 
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Fig. 19 Weighted electronic bands of rutile VO2. The width of the bars given for each band 
indicates the contribution due to the 3d x 2 _ y i orbital of the V atom at (0,0,0) relative to the 
local rotated reference frame. 




Fig. 20 Weighted electronic bands of rutile VO2. The width of the bars given for each band 
indicates the contribution due to the 3d yz orbital of the V atom at (0,0,0) relative to the 
local rotated reference frame. 
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Fig. 21 Weighted electronic bands of rutile VO2. The width of the bars given for each band 
indicates the contribution due to the 3dx Z orbital of the V atom at (0,0,0) relative to the 
local rotated reference frame. 



our work on M0O2 |23| : In all three figures, each band at each k-point is given a bar, 
which length is a measure for the contribution from a specified orbital. Since we again 
refer to the local frame of reference, Figs, [h] to correspond to the partial V 3d ti g 
DOS shown in Fig. which likewise made use of the local coordinate system. 

As in M0O2 and NbC>2 we observe in Fig. [l9| a strong dispersion of the metal 
2>d x i_ y i derived bands along all lines parallel to T-Z. In contrast, the dispersion per- 
pendicular to this line is almost negligible. Again, we attribute this dispersion to the 
overlap of d x i_ y i orbitals at metal sites neighbouring in z direction. As compared to 
M0O2 the width of these bands relative to the total t^g band width is reduced. This 
is due to the increased value of the rutile c axis and, hence, the larger metal-metal 
distance in this direction. In contrast, the smaller value of the a axis in VO2 causes a 
larger interchain interaction. This leads to an increased splitting of the d x 2_ y 2 bands 
along the lines T-M and A-Z as well as along the lines where the strong dispersion is 
observed. 

The V 3d yz derived bands are highlighted in Fig. They show a large dispersion 
in all directions and, as in M0O2, a large splitting at the r-point. The latter is a 
consequence of the underlying body-centered tetragonal lattice formed by the metal 
atoms and results from backfolding the bands from the Brillouin zone of this lattice 
to the smaller one of the simple tetragonal lattice. This effect can be read off from 
the bands along the lines connecting the V- point to the X-, M-, A-, and Z-point, 
respectively. 

Finally, the V Zd xz derived bands display the smallest dispersion. This is due to 
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the fact that these orbitals induce no substantial er-type overlap with neighbouring 
metal atoms. Instead, they mediate the small coupling between the octahedral chains. 
Still, the larger width of these bands as compared to M0O2 results from the increased 
interchain coupling coming with the smaller value of the a lattice constant. 

Of course, the previous analysis of the t2 g bands with its nearly unique assignment 
of a single character to each band as well as the subsequent interpretation in terms 
of the geometric coordination of the atoms and their orbitals benefitted from the 
very small mixing among these states. Exceptions are the hybridization of the d x 2_ y 2 
bands with the d yz - and d^-derived bands along the symmetry line X-R and, to 
a lesser degree, along the line T-A. Nevertheless, we point out that the three ty, g 
bands disperse almost completely independently and are coupled mainly via charge 
conservation. For this reason the dy band may be regarded as a one-dimensional band 
in a three-dimensional embedding background of tt* bands. The situation is thus 
identical to that found in rutile NbC>2 as well as in hypothetical rutile M0O2 |^||, |24| 

Eventually, the finding of an only weak hybridization of the different types of 
bands has a strong impact on our further understanding since it will allow to establish 
a relation between particular atomic displacements and the response of the electronic 
states. This way, we may find a route to an understanding of the mechanisms leading 
to the low-temperature monoclinic structure. 



6 Results for insulating M1-VO2 

6.1 Band structure and density of states 

The electronic bands of monoclinic VO2 in the energy range of the V 3d states are 
displayed in Fig. ^2] along selected high symmetry lines within the simple monoclinic 
Brillouin zone, Fig. |ll](b). The dominant partial densities of states (DOS) are given 
in Fig. ^3]. As before low lying oxygen 2 s states are not included. 

In Figs. ^ and wc identify the same groups of bands as for rutile VO2. They 
comprise 24 oxygen 2p dominated bands well below the Fermi level and two groups 
of 12 and 8 bands, respectively, at and above Ep. These latter two groups trace back 
mainly to V 3d states. Finally, p-d hybridization leads to V 3d and O 2p contributions 
of the order of about 10% in the energy regions where the respective other partner 
dominates. The energetical separation between the O 2p and V 3d derived bands 
turns out to be slightly smaller than in the rutile phase. This is in agreement with 
the XPS results by Blaauw et al. MQ. Again, these general findings conform well with 



the molecular orbital picture drawn in Sec. 5.1. 

Nevertheless, on going into more detail we find distinct differences between the 
results for the rutile and the monoclinic structure. In particular, we witness the 
evolution of a sharp peaks in the V 3d t^ g partial DOS just below Ep and at « 1.3 eV. 
In general, spectral weight is shifted from the central region of this group of bands to 
the edges. As a consequence, the density of states at the Fermi energy is lowered as 
compared to the rutile value and amounts to 2.08 states/f.u./cV. 

In order to understand these differences we proceed displaying in Fig. |24| the single 
symmetry components of the V 3d partial DOS. In doing so we include only the single 
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Fig. 22 Electronic bands of monoclinic VO2 along selected symmetry lines within the simple 
monoclinic Brillouin zone, Fig. F3(b). 
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Fig. 23 Partial densities of states (DOS) of monoclinic VO2 per unit cell. 
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Fig. 24 Partial V 3d t^g and e g densities of states (DOS) of monoclinic VO2. Selection of 
orbitals is relative to the local rotated reference frame. 
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vanadium atom near the corner of the underlying rutilc cell and use the local rotated 
reference frame as defined in Sec. 



3.1. As in Fig. we recognize the almost perfect 



energetical separation of the 3d t^g and e g groups of bands, which results from the 
octahedral crystal field splitting. The small t2g-e g configuration mixing again points 
to the slight distortions of the octahedra. 

Concentrating on the t 2g partial DOS, we find close similarities of the d xz and 
d yz partial DOS to those calculated for the rutile structure given in Fig. [l^. Yet, 
both bands have experienced a considerable energetical upshift and are left almost 
unoccupied. 

While the d xz and d yz partial DOS still resemble each other the d x z_ y 2 partial 
DOS display striking deviations from both these two curves and the results obtained 
for the rutile structure. Comparing Fig. |l2] to Fig. ^ we recognize strong splitting 
of the d x 2_ y 2 partial DOS into two peaks, which are located just below Ep and at 
~ 1.3 cV. In addition, the partial DOS in between is considerably reduced. This leads 
to the decrease of the density of states at the Fermi energy. 

Our findings for the ti g states of monoclinic VO2 are in good agreement with 
the band scheme proposed by Goodenough ||, which we sketched in Fig. g, We 
thus attribute the striking behaviour of the d x 2_ y 2 (or d\\) bands to the pairing of 
the vanadium atoms along the rutile c axis and the resulting bonding-antibonding 
splitting. In contrast, the energetical upshift of the d xz and d yz (or it*) bands results 
from the lateral antifcrroelectric displacement of the vanadium atoms, which increases 
the p-d bonding. 

The situation becomes clearer from a comparison to the band structure of rutile 
VO2 shown in Fig. [2q. Here we have used the first Brillouin zone of the simple 
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Fig. 25 Electronic bands of rutile VO2 along selected symmetry lines within the simple 
monoclinic Brillouin zone, Fig. |lC 
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monoclinic lattice, Fig. g_0|(b) and folded back all bands from the larger Brillouin zone 
of the tetragonal lattice, Fig. |l^(a). Just in the same manner as for M0O2, those 
two bands in Fig. which give rise to the lowest states within the t2g-derived group 
of bands at the T-point, bend upwards along the line T-A and cross the higher lying 
bands. In monoclinic VO2 these two bands are separated from the other states and 
form a split-off doublet, which, however, is still touching the higher lying bands at 
the A-point. For this reason, we end up with a finite density of states at the Fermi 
energy, which spoils the experimental finding of an insulating ground state with an 
optical band gap of about 0.6 eV. Instead, from the touching bands at the A-point and 
the occupied band width of the conduction band at the C-point of about 0.1 eV we 
would rather regard monoclinic VO2 as a semimetal, in agreement with the results by 
Wentzcovitch et al. M . Like these authors we attribute this failure of the calculations 
to the shortcomings of the local density approximation, which usually underestimates 
the band gap by about 50% and in some small gap semiconductors like Ge misses 
the gap at all. This seems to be also the case for monoclinic VO2. Like in the 
zincblende-type semiconductors the optical band gap originates to a large part from 
bonding-antibonding splitting of hybridized bands, which, obviously, is not well enough 
accounted for by the LDA. 

6.2 Chemical bonding 

Again we investigate the chemical bonding by determining the covalence energy. The 
result is shown in Fig. The E cov curves for the O 2p-0 2p and the V 3d-0 
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Fig. 26 Partial covalence energies (E cov ) of monoclinic VO2. 

2p bonding agree fairly well with the respective curves for rutile VO2 as shown in 
Fig. [H| Qualitative differences show up in the V 3d-V 3d bonding. In particular, 
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we observe a change in sign in the energy region of the bonding V 3d x 2_ y 2 bands. 
While the respective E cov curve in the lower half of the V 3d t^ g bands was positive 
(antibonding) for rutilc VO2 it has become negative (bonding) just below the Fermi 
energy in the Mi structure. As for monoclinic M0O2 this signals the increased metal- 
metal bonding due to the pairing of the metal atoms. Furthermore, the formation 
of bonding and antibonding V 3d x 2_ y 2 states below and above the Fermi energy, 
respectively, decreases the value of the total integrated E cov at the Fermi level. 

6.3 Comparison to experiment 

As for metallic VO2 we display in Figs. |27| and total and partial V 3d and O 2p 
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Fig. 27 Total and partial densities of states (DOS) of monoclinic VO2 folded with a 0.25 eV 
wide Gaussian (lower set of curves) and UPS spectra (upper curve; note the offset introduced 
in order to distinguish experimental and theoretical results; from Ref. | p3| , ptHf). 

densities of states. The curves are folded with a Gaussian of 0.25 and 0.5 eV width, 
respectively, for the occupied and unoccupied part of the spectrum. Good agreement 
with the XPS and UPS spectra by Blaauw et al., Shin et al., Bermudez et al., and 
Goering et al. is found ||, EH HI II H- Data h Y the latter 

group were added in 

Fig. ^?]. According to these experiments the UPS spectra taken in the metallic and 
the insulating phase are very similar with respect to the oxygen dominated bands but 
show distinct differences in the low-binding region. They arise from a sharpening and 
slight downshift of the V 3c?-derived peak |4l], |4^, As a comparison of Fig. |2?] 
with Fig. [l4| reveals, this experimental result is well accounted for by our calculations 
except for the fact that the V 3d downshift is underestimated. 

More striking differences between metallic and insulating VO2 are observed in the 
soft-X-ray absorption spectra (XAS) as measured by Miiller et al. ^8|, which we 
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Fig. 28 Partial O 2p densities of states (DOS) of monoclinic VO2 folded with a 0.5 eV wide 
Gaussian (lower set of curves) and XAS O K edge spectra (upper set of curves; note the offset 
introduced in order to distinguish experimental and theoretical results; from Ref. |)^, ^8j, 
data shifted by 529.5 eV). 



display in Fig. [28|. The experimental data were shifted by 529.5 eV. As already de- 
scribed in Sec. 5A, the XAS spectra allowed to probe the angular dependence by 
varying the orientation of the polarization vector E from parallel (</> = 90°) to perpen- 
dicular (<fi = 0°) to the rutile c axis. Dipole selection rules then allow for transitions 
from O Is to cither the O 2p z or the O 2p x and 2p y states. This is reflected by the 
calculated partial DOS of the corresponding final states. However, note that the curve 
marked O 2p x , p y actually comprises only the mean average of these orbitals since the 
component perpendicular to E is not seen in experiment. For the same reason, the 
curve marked O 2p actually contains the full contribution from the p z orbital but only 
half of the contributions from the p x and p y orbitals. 

According to the experiments by Miiller et al. as well as by Abbate et al. major 
differences between the XAS spectra of metallic and insulating VO2 originate from 
the appearance of an additional peak in the low- temperature phase |9(], ^|. In 
the spectra by Abbate et al. this new peak is located at 530.9 eV and seems to have 
split off from the peak at 529.9 eV, which, like the one at 532.5 eV, exists also in the 
metallic phase jl7|. Abbate et al. interpreted the additional peak as arising from an 
upshift of the antibonding d\\ band by about 1 eV. The other two t2 g bands, namely 
the 7T* bands, stay at their original position. The XAS curves taken by Miiller et al. 
are almost identical to the results of Abbate et al.. However, in a more detailed fit 
Miiller assigned the spectra to five different peaks at 529.9, 530.4, 531.2, 532.7, and 
533.9 eV @. 

Again, there is an overall good agreement between experiment and theory with 
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respect to positions, intensities, and angular dependencies of the peaks. As a compar- 
ison with Fig. [[5] reveals, the calculated curves reflect the upshift of the O 2p z peak 
from 0.7 cV to 1.3 eV, which is the position of the upper d x 2_ y 2 peak in Fig. At 
the same time, all other peaks hardly move. Nevertheless, the O 2p z peak has its 
maximum at the same energy as the O 2p x ,p y curve at 1.3 eV. This is in contrast to 
the XAS experiments, which find the maximum of the O 2p z peak near 1.8 eV. Thus, 
it seems that the upshift of the O 2p z states and, hence, of the upper d x 2_ y 2 states, 
is underestimated by the calculation in the same manner as in our previous work on 



M0O2 1 23 . From a comparison of the calculated DOS with both angle-resolved UPS 
and XAS data for the molybdenum dioxide, we obtained perfect agreement with ex- 
periment except for the bonding and antibonding Mo 4d x 2_ y 2 bands. Their splitting 
was underestimated by about 1 eV. For VO2 the calculations yield a separation of 
w 1.7 cV between the bonding and antibonding d x 2_ y 2 bands. A value of 2.0 to 2.5 eV 
has been extracted from UPS and XAS measurements 41, 47|. 

Finally, as for rutile VO2 the calculated high-energy peak at sw 4.8 eV gives rise to 
only a very slight shoulder in the experimental curve. Yet, it can be fully resolved by 
the aforementioned five peak analysis by Muller p3. 



6.4 Embedded Pcicrls instability 

The upshift of the V 3d xz and d yz bands as well as the bonding-antibonding splitting 
of the d x 2_ y 2 band as observed in the partial DOS already signaled the general validity 
of the mechanism proposed by Goodenough and, more recently, by Wcntzcovitch et 
dl. [[|, 0, H . Nevertheless, we still aim at obtaining more direct support for these ideas 
from the electronic structure. In order to prepare for the discussion we display in Fig. 
p9| the weighted band structure of rutile VO2 corresponding to the d x 2_ y 2 states. It 
results from that shown in Fig. [l^ by backfolding to the monoclinic Brillouin zone, 
Fig. |l^(b) . As before the length of a bar given with a band at a k-point is a measure 
of the contribution from the d x 2_ y 2 orbital. Note that again we refer to the local 
reference frame. 



In Fig. 29 we observe an arrangement of the d x 2_ y 2 bands in two parts below and 
above 0.4 eV, respectively. They are connected by dispersing bands especially along 
the lines T-A and T-B. Both lines have components perpendicular to the rutile basal 
plane and correspond to half the line T-R of the tetragonal Brillouin zone. In contrast, 
the dispersion parallel to the rutile basal plane is suppressed to a large degree. The 
respective bands stay at energies either well below or above 0.4 eV, even if there is some 
dispersion within the lower group. Finally, d x 2_ y 2 bands along the lines A-E and B- 
D, while likewise displaying a reduced dispersion, stay at an intermediate energy in 
the middle between the two subgroups. Note that these lines lie completely within 
horizontal planes at height ^ . The general characteristics of these bands are thus very 
similar to those of the Ad x 2_ y z bands of hypothetical rutile M0O2. In that system the 
dispersionless bands along A-E and B-D were responsible for the flat portions of the 
Fermi surface. They were thus interpreted as driving the Peierls-type instability of 
this structure. 

The V 3d xz and d yz bands can be identified in Fig. ^ as the states without a bar. 
This interpretation is based on the fact that the O 2p and the V 3d e g states play 
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Fig. 29 Weighted electronic bands of rutile VO2. The width of the bars given for each band 
indicates the contribution due to the 3d x 2 _ y i orbital of the V atom at the corner of the rutile 
cell relative to the local rotated reference frame. 



only a minor role in the energy range shown. In contrast to the d x 2_ y 2 bands these 
7r* bands do not fall into energetically lower and upper bands. Instead they rather 
show isotropic dispersion. Finally, we note again the much reduced hybridization 
between both types of bands, which we already mentioned at the end of Sec. 5.6. As a 
consequence, the d|| and n* bands can, as before, be regarded as nearly independent. 

The corresponding weighted band structure for monoclinic VO2 is shown in Fig. ||(]. 
The differences between both figures is striking. On going from rutile to monoclinic 
VO2 we observe a strong splitting of the d x 2_ y 2 bands into two narrow subbands at 
the lower and upper edge of the t^g group of bands. This splitting is a consequence of 
the metal-metal dimerization. In the partial DOS, Fig. ^J, this bonding-antibonding 
splitting led to the formation of the two distinct peaks just below Ef and at ps 1.3 eV. 
In particular, the d x 2_ y 2 bands along the lines T-A-E-Z and T-B-D-Z have moved 
downward and upward, respectively, by about 0.5 eV relative to their positions in 
the rutile structure. As a consequence, the so far unoccupied branch of these bands 
between Ep and « 0.4 cV has become almost fully occupied and the Fermi surface 
due to the d x 2_ y 2 states has almost vanished. Moreover, these changes explain the 
occurrence of the split-off double band at the lower edge of the ti g group of bands, 



which we mentioned at the end of Sec. 6.1. 

In contrast, the V 3d xz and d yz bands, i.e. those bands in Fig. 30, which have 



bars of negligible lengths, experience a considerable upshift by fts 0.5 eV due to the 
zigzag-like antiferroelectric displacement of the vanadium atoms. This result is in full 
agreement with the data by Shin et al. Ell . Yet, with only minor exceptions, the shift 
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Fig. 30 Weighted electronic bands of monoclinic VO2. The width of the bars given for each 
band indicates the contribution due to the 3d x 2 orbital of the V atom at the corner of the 
rutile cell relative to the local rotated reference frame. 



of the ir* bands leaves the general topology of these bands intact and causes their 
almost complete depopulation. Our findings agree well with the arguments given by 
Goodenough, who concluded from the isotropic conductivity that the ir* states should 
be located at the lower edge of the conduction band well below the antibonding d\\ 
band R]. Furthermore, Goodenough deduced from the measured low Hall mobility 
a narrowing of the it* bands |Q. Indeed, our calculations result in a band width of 
about 2 eV of the it* bands, which is considerably less than the 2.5 eV found in the 
rutile phase. 

Finally, we note that the hybridization between both types of bands is still very 
small. In other words, the just described changes of the d\\ and ir* bands occurred 
rather independently. The coupling between them is still only via the common Fermi 
energy. Just in the same manner as in Nb02 and M0O2 we may thus regard the 
transition to the low-temperature state as a Peierls instability of the d|| bands in an 
embedding background of tt* electrons. 

An aspect, which seemingly has not been considered before, arises from the fact, 
that metal-metal pairing along the rutile c axis, in addition to increasing the bonding- 
antibonding splitting of the d x 2_ y 2 bands, shifts the vanadium atoms also relative 
to the surrounding oxygen octahedra and, hence, causes an increased metal-oxygen 
overlap. Thus, if not compensated by a corresponding displacement of the oxygen 
octahedron, the dimcrization likewise brings in an antiferroelectric component. As a 
consequence, within an investigation of the monoclinic Mi structure we are actually 
not able to uniquely attribute the upshift of the tt* bands to either the zigzag-type in- 
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plane displacement or the metal-metal pairing. This is due to the fact that in the Mi 
phase all metal chains experience the same atomic displacements. The issue can only 
be resolved by studying the monoclinic M2 phase, which brings in a disproportionation 
of the vanadium chains. 



7 Results for insulating M2-VO2 

In the previous two sections we have investigated in detail both the metallic rutile 
and the insulating monoclinic (Mi) phase of V02- In addition, we proposed a possible 
explanation for the origin of the metal-insulator transition in terms of an embedded 
Peierls-like instability. Our conclusions confirmed the previous findings of Wcntzcov- 
itch et al., who, using a molecular dynamics type approach, arrived at very similar 
results |8| pi. Still, the discussion by these authors was criticized by Rice et al, who 
focused attention on the monoclinic M2 phase of VO2 Jl3| . 

The arguments given by Rice et al. aimed at attempts to explain the metal-insulator 
transition of stoichiometric VO2 purely in terms of the dimcrization in the monoclinic 
Mi phase. If this were the only or the dominating driving force there would be no 
reason why VO2 in the M2 phase would be insulating, since only half of the chains 
dimcrize. Yet, as has been already pointed out by Goodenough, metal-metal pairing 
is not the only source of an efficient splitting of the d\\ band but the same effect could 
be achieved by exchange splitting coming with the antiferromagnetic order @, . The 
latter was indeed observed on the zigzag chains of the M2 structure. Nevertheless, 
from the antiferromagnetic ordering several authors concluded that the respective 3d 
electrons would be rather localized and, while being accessible to strong electron- 
electron correlations, resist a proper description within band theory. Due to the close 
relation between the low-temperature monoclinic phases these arguments would then 
also hold for the Mi phase. 

Since the M2 phase has not yet been considered using state-of-the-art electronic 
structure calculations we included its investigation in the present study. From discus- 
sion of the M2 phase we expect a deeper understanding of the mechanisms, which are 
active in the Mi phase, since each chain shows only one type of octahedral deforma- 
tion. Furthermore, we should be able to distinguish the different types of vanadium 
atom displacements with respect to their coupling to the electronic states and, in par- 
ticular, to identify the implications of the two different types of antiferroelectric modes 
arising from the zigzag-type and pairing displacements. The latter can not be resolved 
in investigations of the Mi phase. 



As has been outlined in Sec. 2^3, Pouget et al. proposed to view the monoclinic M2 
phase as a metastable modification of the Mi phase and to regard both the doping 
with Cr or Al and the application of uniaxial stress as only small perturbations of 
pure VO2 under ambient pressure. These arguments were supported by the crystal 
structure data reported for different types of dopants. For both the rutile and the M2 
phase the parameters given by Marezio et al. for Vo.976Cro.o2402 and by Ghedira et 
al. for V0.985AI0.015O2 are almost identical p0| , 66 and for the rutile structure they 



arc identical to the numbers given for stoichiometric VO2, see Sec. B.l and Tabs, 
and ^. It is thus well justified to use the crystal structure data for the doped material 
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as input for the present study of the M2 phase. 

Our investigations for the M2 phase proceeded in two steps. First we performed 
calculations using the crystal structure data given by Marezio at at and assuming 
spin-degeneracy. Spin-polarization was included in a second step. This way we were 
able to check the stability of the antifcrromagnctic state with respect to the spin- 
dcgcncrate case. Furthermore, our procedure allowed to separate the impact of the 
structural transformations and magnetic order on the electronic structure. 




Fig. 31 Electronic bands of monoclinic M2 VO2 along selected symmetry lines within the 
simple monoclinic Brillouin zone. 

the centered monoclinic Brillouin zone. Except for the small lattice strain, which 
distorts the basal plane, the latter Brillouin zone is identical to the simple monoclinic 
Brillouin zone of the Mi phase. The dominant partial densities of states (DOS) are 
given in Fig. [52|. As for the Mi phase of VO2 we observe many similarities with the 
rutile results as concerns the ordering, widths and separations of the bands. Still, we 
find the O 2p-, V 3d t% g - and e 9 -dominated groups of bands below, at and above the 
Fermi energy, respectively, as well as a hybridization of these states of the order of 
10%. 

Differences between the results for the two low-temperature phases show up on 
closer inspection of the V 3c? derived partial DOS. Due to the presence of two different 
chains in the M2 phase we have distinguished the two types of vanadium atoms, namely 
Vi and V2, which belong to the dimerizing and the zigzag chains, respectively. Distinct 
differences occur in the t2 g group of bands. The Vi 3d ti g partial DOS displays strong 
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Fig. 32 Partial densities of states (DOS) of monoclinic M2 VO2 per unit cell. 



splitting into two peaks just below Ep and at w 1.8 cV similar to that observed for 
the Mi phase, Fig. ||^. In contrast, the V2 3d ti g partial DOS is dominated by a 
single maximum at « 0.3 eV and resembles the respective partial DOS found for the 
rutile structure, Fig. fll Qualitatively, the same holds for the V 3d e g partial DOS, 
for which we note the close similarity of the partial DOS assigned to the Vi and V2 
atoms to those calculated for the Mi and the rutile structure. All these findings reflect 
our expectations. In a non-spinpolarized calculation the dimerizing and zigzag chains 
should display similar behaviour as that known from the Mi and the rutile phase, 
respectively. 

Going into more detail we display in Figs. |3^ and |34| the 3d partial DOS of the 
Vi and V2 atoms as decomposed into their symmetry components relative to the 
respective local reference frame of each chain. As for the rutile and the Mi phase we 
observe for both types of vanadium atoms the energetical separation of the 3d ti g and 
e g groups of bands resulting from the crystal field splitting. However, within these 
groups of bands we also witness the aforementioned differences between the Vi and 
V2 atoms. 

The t2 g partial DOS of the paired Vi atoms closely resemble the ti g partial DOS 
of the Mi phase as given by Fig. |4| They both are characterized by the strong 
bonding-antibonding splitting of the d x 2_ y 2 states, which is even more pronounced in 
the M2 phase. This result fully agrees with the NMR data, which revealed the strong 
similarity of the Knight shift due to the dimcrized vanadium atoms to that measured 
in the Mi phase (6^]. In contrast, the d x i_ y i partial DOS of the V2 atoms, which 
belong to the zigzag chains, display a completely different behaviour. These states do 
not show any splitting but form a single maximum at w 0.4 eV with a broad shoulder 
extending to -0.65 cV and only minor contributions at energies above w 0.7. This 
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Fig. 33 Partial Vi 3d ti g and e g densities of states (DOS) of monoclinic M2 VO2. Selection 
of orbitals is relative to the local rotated reference frame. 



shape thus bears a strong similarity to the corresponding DOS of the rutile phase, 
see Fig. 12, rather than to that of the monoclinic Mi phase. To conclude, the d x i_ y i 
partial DOS can be readily understood from the local environment of the respective 
atoms as being more Mi- and rutile-like, respectively. 

The d xz and d yz curves follow this general trend in that the Vi partial DOS reflect 
the behaviour known from the Mi phase while the V2 partial DOS are more rutile- 
like. Yet, two features are worth mentioning: i) The splitting of the 3d yz partial 
DOS as observed in the rutile structure is almost completely suppressed on the V2 
chains. This is due to the zigzag-like displacement of these atoms, which undermines 
the in-plane overlap of the d yz orbitals and, hence, reduces their bonding-antibonding 
splitting, ii) The upshift of the ir* states affects both the Vi and the V2 partial DOS. 
From the discussion in Sec. |[ where we attributed the upshift of these bands to the 
antifcrroclcctric mode of the vanadium atoms coming with the zigzag-like in-planc 
displacements, we would have expected a raise in energy only for the ir* states of the 
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Fig. 34 Partial V2 3d ti g and e 9 densities of states (DOS) of monoclinic M2 VO2. Selection 
of orbitals is relative to the local rotated reference frame. 



V2 atoms but not for the vanadium atoms in the dimcrizing chains. Yet, a closer look 
at the partial DOS in Figs. [33] and [m| reveals that the centers of gravity of the 7r* 
bands actually is higher on the dimerizing chains. 

An explanation for this latter behaviour can be found in an analysis of the possible 
metal-oxygen overlaps using Tab. where we combine those particular V 3d ti g and 
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O 2p orbitals. which, within a molecular-orbital picture, are expected to overlap. We 
complement the table with the partial O 2p partial DOS in Fig. [35], where we have 
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Fig. 35 Partial O 2p densities of states (DOS) of monoclinic M2 VO2. Selection of orbitals 
is relative to the local rotated reference frame of atom Vi. 
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again used the local rotated frame of reference for the selection of the orbitals. Note 
that we opted for the local reference frame corresponding to the dimcrizing Vi chains. 



As outlined in Sec. 3.1 interpretation in terms of the reference frame of the zigzag-like 
V2 chains simply requires exchange of the y and z components. For the following 
discussion we recall that the Oi atoms are the apical oxygen atoms of the dimcrizing 
Vi chains and occupy the equatorial positions of the zigzag-like V2 chains, see Fig. [?]. 
For oxygen atoms O2 and O3 things are reversed with the latter having the shorter 
distance to the central V2 atom. In addition, due to the dimerization, O2 have a very 
short bond with the vanadium atoms of type Vi . 

The aforementioned geometry is reflected by the partial DOS. In Fig. [35] we observe 
a striking similarity of the partial DOS corresponding to atom O2 to the Vi 3d partial 



DOS as shown in Fig. 33. In general, all contributions in Figs. |33| to |35| can be well 
understood using the orbital overlaps listed in Tab. |l3|. In particular, we observe the 
hybridization of the V2 3d xz peak at ss 0.4 eV with the O3 2p x contribution. Onset of 
the zigzag-like mode increases this overlap and shifts these states to higher energies. 
In contrast, the Vi 3d xz states, which are found mainly above 0.8 eV, overlap with 
the 2p z states of their equatorial O2 and O3 partners. This overlap increases with 
increasing dimerization on these chains and, hence, also causes an upshift in energy. 
We thus witness a considerable antiferroclcctric component of the metal-metal pairing, 
which drives upshift of the 7r * st ates on the corresponding chains and which we already 



proposed at the end of Sec. |6.4 . 

The interpretation of the results found especially for the d x i_ y i partial DOS of the 
two different types of chains is confirmed by the weighted band structures as displayed 
111 Figs. |6| and |7|. A gain we note the energetical upshift of the d xz and d yz bands, 
which is almost independent of the type of chain, on which these orbitals are located. 
The d x 2_ y 2 bands of the Vi atoms show the strong bonding-antibonding splitting 
already observed in the partial DOS, Fig. in complete analogy to the weighted 
d x 2_ y 2 bands of the Mi phase, Fig. |3(]. In both phases it results from the metal-metal 
pairing on the respective chains and leads to the semimctallic behaviour of these bands 
with a clear separation of bands of different characters. The same effect was already 
observed in the results for the Mi structure. In contrast, the d x 2_ y 2 bands of the V2 
atoms, which have equal metal-metal distances, resemble the dispersion of these bands 
in the rutile structure, see Fig. ^9|. Their reduced band width, which is by « 0.4 eV 
smaller as compared to the rutile structure, reflects the increased V2-V2 distance in the 
zigzag chains. Nevertheless, these bands still display the one-dimensional dispersion 
found for the rutile phase. As a consequence, from the spin-degenerate calculations we 
obtain half-filled V2 d x 2_ y 2 bands of w 1.1 cV width, which are clearly metallic. To 
conclude, enforcing spin-degeneracy we find indications for a metal-insulator transition 
very similar to those discussed in Sec. || in the dimerizing chains but at the same time 
we obtain rutilc-likc metallic behaviour due to the bands originating from states on 
the zigzag chains. 



7.2 Spinpolarizcd calculations 



Next we turn to the spinpolarized calculations as suggested by the experimental finding 
of antiferromagnetic Heisenbcrg chains in the M2 phase of VO2. As a result we obtain 
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Fig. 36 Weighted electronic bands of monoclinic M2 VO2. The width of the bars given for 
each band indicates the contribution due to the 3d x 2 _ y i orbital of atom Vi, which belongs 
to a dimerizing chain, relative to the local rotated reference frame. 
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Fig. 37 Weighted electronic bands of monoclinic M2 VO2. The width of the bars given for 
each band indicates the contribution due to the ?>d x i_ y i orbital of atom V2, which belongs 
to a zigzag chain, relative to the local rotated reference frame. 
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a stable antiferromagnetic solution with a total energy, which is by 0.7mRyd per 
unit cell lower than that growing out of the non-spinpolarized calculation. Magnetic 
moments of 0.46 fiB are carried by the V2 3d states of the zigzag chains. The moments 
of all other atoms arc cither negligibly small or, for symmetry reasons, exactly zero 
as in case of the Vi. The dominant partial DOS are displayed in Fig. pq. Note 
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Fig. 38 Partial densities of states (DOS) of antiferromagnetic monoclinic M2 VO2 per unit 
cell. 

that all curves in Fig. |38| comprise contributions from one magnetic sublattice only. 
On comparing Fig. [58] to Fig. [52] for the non-spinpolarized case we recognize the main 
changes arising from spinpolarization in the small energy interval from -1 to 1 eV about 
the Fermi energy. The Vi 3d t 2g partial DOS are identical for both spin directions due 
to the vanishing magnetic moments at these sites. Spin splitting arises mainly from 
the V2 3c? t2 g states. Their spin up and down DOS arc dominated by a broader peak 
just above Ep and rather sharp peak at ~ 0.7 eV, respectively. 

The changes originating from antiferromagnetic order are reflected by analysis of 
the V 3d partial DOS, which we display in Figs. [59] and [h)] Again, the identical spin 
up and down DOS of the Vi atoms reflect the vanishing magnetic moments of these 
atoms. In contrast, we observe spin-dependent shifts of the V2 partial DOS in the 
energy interval between -1 and 1 eV and, to a lesser degree, between 2 and 4eV. While 
spin-splitting of the d xz and d yz states is rather small on these atoms the magnetic 
moment originates mainly from the local d x 2_ y 2 states. 



Recalling the non-spinpolarized d x 2_ y 2 partial DOS as shown in Fig. 34 with its 
prominent peak at w 0.5 eV and the broad shoulder below Ep we can now attribute 
the peak at 0.5 eV mainly to the spin down states while the broad shoulder comprises 
almost exclusively the spin up states. This situation becomes clearer from the weighted 
electronic structures as shown in Figs, fo and E4 These band structures highlight 
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Fig. 39 Partial Vi 3d t2 9 and e 3 densities of states (DOS) of antiferromagnetic monoclinic 
M2 VO2. Selection of orbitals is relative to the local rotated reference frame. 



the spin up and down V2 3d x 2_ y 2 bands relative to the local reference frame of the 
single atom belonging to the spin up sublattice. While the occupied 3d x 2_ y 2 bands 
are mainly of spin up character, bands at higher energies are almost exclusively of 
spin down character. The spin majority and minority bands thus play a similar role 
as the bonding and antibonding bands of the vanadium 3d x 2_ y 2 states located at the 
dimcrizing chains. The spin majority states arc almost completely filled and show only 
small contributions to the bands at w 0.6 cV. In contrast, the spin minority states 
dominate the bands in the energy range between 0.4 and 0.6 eV. As a consequence, 
antiferromagnetic ordering is obtained in the zigzag chains. At the same time, singlet 
formation arises in the dimerizing chains. 

As for the Mi phase, the band structures and densities of states calculated with 
the M2 structure lack the observed optical band gap. Yet, the situation is very similar 
to that found in Sec. EL We find an overlap of the valence and conduction bands of 
w 0.2 eV, which we attribute again to the shortcomings of the LDA. However, bands 
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Fig. 40 Partial V2 3d £23 and e g densities of states (DOS) of antiferromagnetic monoclinic 
M2 VO2. Selection of orbitals is relative to the local rotated reference frame. 



of different orbital character are clearly separated. For this reason, it is well justified 
to regard DFT-based calculations as being capable to describe also the metal-insulator 
rutile-to-M2 transition. This holds even if the opening of the optical band gap is just 
missed due to the LDA. 

In summary, while the rutile phase is characterized by two types of rather weakly 
hybridized quasi-oncdimensional dii and isotropically dispersing n* bands the changes 
in crystal structure and magnetic order coming with the M2 structure lead to dis- 
tinct changes in the electronic structure. Both the zigzag-type displacement on the 
V2 chains and the dimerization on the Vi chains induce antifcrroclectric shifts of the 
vanadium atoms inside the surrounding oxygen octahedra. The resulting increase of 
metal-oxygen bonding shifts the ir* bands to higher energies. In addition, dimer- 
ization of the Vi atoms causes splitting of the corresponding dn bands into bonding 
and antibonding branches. In contrast, the d\\ bands belonging to the V2 chains dis- 
play spin-splitting and antiferromagnetic order in these chains. As for the Mi phase 
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Fig. 41 Weighted electronic bands of antiferromagnetic monoclinic M2 VO2. The width of 
the bars given for each band indicates the contribution due to the spin up 3d x 2_ y 2 state of 
atom V2, which belongs to a zigzag chain, relative to the local rotated reference frame. 
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Fig. 42 Weighted electronic bands of antiferromagnetic monoclinic M2 VO2. The width of 
the bars given for each band indicates the contribution due to the spin down 3d x 2_ y 2 state 
of atom V2, which belongs to a zigzag chain, relative to the local rotated reference frame. 
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coupling between the different types of bands is still weak and mainly via the com- 
mon Fermi energy. The scenario of an embedded Pcicrls instability proposed for the 
transition from the rutile to the Mi phase has to be extended by the possibility of 
antifcrromagnctic ordering as already mentioned by Goodenough. 



8 Conclusion 

In the present first principles study, augmented spherical wave calculations were used 
to investigate the structural, electronic, and magnetic properties as well as the metal- 
insulator transitions of vanadium dioxide. The essential features of the metallic rutile 
as well as the insulating Mi and M2 phases are well described within DFT and LDA. 
In metallic VO2, the electronic structure consists of oxygen 2p bands well below the 
Fermi energy, as well as crystal field split vanadium 3c? ti g and e a g bands. States near 
Ep consist of two very weakly hybridizing types of ti g bands, namely, the d\\ states with 
a quasi-onedimensional dispersion parallel to the characteristic V chains, and isotropi- 
cally dispersing e g bands. In the low-temperature Mi structure of stoichiometric VO2, 
metal-metal dimcrization splits the dn band into bonding and antibonding branches, 
whereas the e g states shift to higher energies due to reduced V-0 distances. Although, 
due to the imperfections of the LDA, semimctallic behaviour with a band overlap of 
« 0.1 cV is obtained rather than the observed optical band gap, the d\\ and e g bands 
are completely separated. Since both types of bands are still coupled by charge con- 
servation rather than hybridization, the insulating state is interpreted as due to a 
Peierls-like instability of the d\\ bands in an embedding reservoir of e g electrons. The 
situation is thus identical to the one previously found for M0O2 and NbC>2 |2||, pj) , 
and gives rise to a unified picture for the transition-metal dioxides at the beginning of 
the d series. This approach explains both (i) the destabilization of the rutile structure 
in terms of increased metal-metal bonding and (ii) the metal-insulator transitions of 
the d 1 members. 

In the insulating M2 phase, obtained on doping or applying uniaxial pressure, which 
generally is viewed as a metastable phase of stoichiometric VO2, the metal- insulator 
transition results from different, complementary mechanisms. Still, depopulation of 
the e g states is due to the antiferroelectric components of both the zigzag-like displace- 
ments and the dimcrization, which increase the metal-oxygen bonding. In addition, 
the embedded Peierls-like instability proposed for the rutile-to-Mi transition causes 
bonding-antibonding splitting of the dn bands connected with the dimerizing chains. 
In contrast, the d|| states located on the chains with equidistant vanadium atoms ex- 
perience spin-splitting due to antifcrromagnctic order of the localized moments, which 
are carried predominantly by the d\\ electrons. As for the other phases, hybridization 
between the and e g states is very weak, and the transition can be interpreted as a 
combined embedded Peierls-like/antiferromagnetic instability of rutile VO2. 
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